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(N ; ABSTRACT 



Context. The iron Ka line at 6.4 keV provides a valuable spectral diagnostic in several fields of X-ray astronomy. The line often 
results from the reprocessing of external hard X-rays by a neutral or low-ionized medium, but it can also be excited by impacts of 
low-energy cosmic rays. 

Aims. This paper aims to provide signatures allowing identification of radiation from low-energy cosmic rays in X-ray spectra showing 
the 6.4 keV Fe Ka line. 

Methods. We study in detail the production of nonthermal line and continuum X-rays by interaction of accelerated electrons and ions 
with a neutral ambient gas. Corresponding models are then applied to XMM-Newton observations of the X-ray emission emanating 
from the Arches cluster region near the Galactic center. 

Results. Bright 6.4 keV Fe line structures are observed around the Arches cluster. This emission is very likely produced by cosmic 
i-C 1 rays. We find that it can result from the bombardment of molecular gas by energetic ions, but probably not by accelerated electrons. 

Using a model of X-ray production by cosmic-ray ions, we obtain a best-fit metallicity of the ambient medium of 1.7 ± 0.2 times 
the solar metallicity. A large flux of low-energy cosmic ray ions could be produced in the ongoing supersonic collision between the 
star cluster and an adjacent molecular cloud. We find that a particle acceleration efficiency in the resulting shock system of a few 
+jf 1 percent would give enough power in the cosmic rays to explain the luminosity of the nonthermal X-ray emission. Depending on the 

unknown shape of the kinetic energy distribution of the fast ions above ~ 1 GeV nucleon -1 , the Arches cluster region may be a source 
of high-energy y-rays detectable with the Fermi Gamma-ray Space Telescope. 

Conclusions. At present, the X-ray emission prominent in the 6.4 keV Fe line emanating from the Arches cluster region probably 
offers the best available signature for a source of low-energy hadronic cosmic rays in the Galaxy. 
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\ 1. Introduction Vink et al. 1 1997ft , and molecular clouds in the Galactic center 

O'tk u v r fHAVxrt t it i - HPt region (Pontietal. [2011. 
The Fe Ka line at 6.4 keV from neutral to low-ionized Fe atoms 

rl ■ is an important probe of high-energy phenomena in various as- ° ne of the best studied cases of this emission is the Sun. The 

; trophysical sites. It is produced by removing a K-shell electron, observed line intensity and light curve during several flares sug- 

. . . either by hard X-ray photoionization or by the collisional ion- § est that excitation of Fe atoms occurs mainly by photoiomza- 

> ■ ization induced by accelerated particles, rapidly followed by an tlon lnduced b y flare X " ra y s wlth ' however, an additional contn- 

electronic transition from the L shell to fill the vacancy. The Fe butlon ln some impulsive events fromcolhsional ionization by 

Ka line emitted from a hot, thermally-ionized plasma at ioniza- accelerated electrons (Zarro et al. [1993. A contribution from 

c3 1 tion equilibrium is generally in the range 6.6-6.7 keV depending collisional ionization by accelerated electrons is also discussed 

on the plasma temperature ^ or tbe emission at 6.4 keV from low-mass flaring stars (Osten 

The 6.4 keV Fe Ka line is a ubiquitous emission feature et "LEHB and young stellar objects (Giardino et al.[2007}. 

in the X-ray spectra of active galactic nuclei (Fukazawa et al. The 6.4 keV line emission from t he Gal actic center (GC) 

1201 11 1. It is also comm only d etected from high-mass X-ray bi- region was predicted by S unyaev et al. dl993l l before being dis- 

naries (Torrejon et al. 1201 Oi l and some cataclysmic variables covered by Koyama et al. dl996b . These authors suggest that the 

(Hellier & Mukai 2004). In these objects, the line is attributed neutral Fe Ka line can be produced in molecular clouds, together 

to the fluorescence from photoionized matter in the vic inity of a with nonthermal X-ray continuum radiation, as a result of repro- 

compact, bright X-ray source (George & Fabian T99T1 and refer- cessed emission of a powerful X-ray flare from the supermassive 

ences therein). The 6 .4 keV line is also detected in sola r flare s black hole Sgr A*. Recent observations of a temporal variation 

(Culhane et al. 1198 11 . low-mass, flaring stars (Osten et al. 1201 Oi l, in the line emission from various clouds of the central molecu- 

massive stars (77 Car; H amagu chi et al. 120071 1. young stellar ob- lar zone can indeed be explained by a long-duration flaring ac- 

jects (Tsujimoto et al. 120051 1. supernova remnants (RCW 86; tivity of Sgr A* that e nded about 1 00 ye ars ago (Mun o et al . 

120071 Inui et al. |2009l Ponti et al. l20T0l Terrier et al. l20T0l l. 

Send offprint requests to: V. Tatischeff Some data also suggest there is a background, stationary emis- 
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sion in the Fe line at 6.4 keV (Ponti et al. 2010), which might 
be due to the interaction of cosmic rays with molecular clouds. 
Observations showing a spatial correlation between the X-ray 
line emission and nonthermal radio filaments have been inter- 
preted as evidence of a large population of accelerated electrons 
in the GC reg ion (Yusef-Z adeh et al. l2002all2"00"7j l. Alternatively, 
Dogiel et al. (2009, 201 1 ) suggest that the neutral or low ioniza- 
tion Fe Ka line from this region could be partly excited by sub- 
relativistic protons generated by star accretion onto the central 
supermassive black hole. 

Low-energy cosmic -ray electrons propagating in the inter- 
stellar medium (ISM) have also been invoked to explain the pres- 
ence of a nonthermal continuum and a weak line at 6.4 keV in the 
spectrum of the Galactic ridge X-ray background (Valina et al. 
2000). But most of Galactic ridge X-ray emission has now been 
resolved into discrete sources, probably cataclysmic variables 
and coronally active stars (Revnivtsev et al. 120091 ). Therefore, 
it is likely that the 6.4 keV line in the Galactic ridge spectrum is 
produced in these sources and not in the ISM. 

In this paper, we study in detail the production of nonthermal 
line and continuum X-rays by interaction of accelerated elec- 
trons, protons, and a-particles with a neutral ambient gas. Our 
first aim is to search for spectral signatures that allow identifica- 
tion of cosmic -ray-induced X-ray emission. We then apply the 
developed models to the X-ray emission from the Arches cluster 
region near the GC. 

The Arches cluster is an extraordinary massive and dense 
cluster of young stars, with possibly 160 O-type stars with ini- 
tial masses greater than 20 M Q and an average mass density 
of ~3 x 10 5 M Q pc~ 3 (Figer et al. [2002b . The X-ray emission 
from the cluster is a mix of thermal and nonthermal radiations. 
The thermal emission is thought to arise from multiple colli- 
sions between strong winds from massive stars (Yusef-Zadeh et 
al. 12002b; Wang et al. 2006). This interpretation was recently 
reinforced by the detection with the XMM-Newton observatory 
of X-ray flaring activity within the cluster, which likely origi- 
nates in one or more extreme colliding wind massive star bina- 
ries (Capelli et al. 12011a]) . Diffuse nonthermal emission promi- 
nent in the Fe Ka 6.4-keV line has also been detected from a 
broad region around the cluster (Wang et al. 2006; Tsujimoto et 
al. 120571 Capelli et al. 1201 lbb . Wang et al. (2006) suggest from 
a 100-ks Chandra observation that this component may be pro- 
duced by interaction of low-energy cosmic -ray electrons with a 
dense gas in a bow shock resulting from the supersonic collision 
of the star cluster with a molecular cloud. In this scenario, the 
nonthermal electrons may be accelerated in the bow-shock sys- 
tem itself and/or in shocked stellar winds within the Arches clus- 
ter. The latter assumption is supported by the detection with the 
Very Large Array (VLA) of diffuse nonthermal radio continuum 
emission from the cluster (Yusef-Zadeh et al. 2003). However, 
Tsujimoto et al. (120071 ) show from Suzaku observations using 
preliminary calculations of Tatischeff (2003) that this scenario 
would require a very high Fe abundance in the ambient medium, 
about four to five times the solar value. Capelli et al. ((201 lb) 
have recently favored a photoionization origin for the 6.4 keV 
line from the Arches cluster region, although not excluding a 
production by low-energy cosmic -ray electrons and/or protons. 
We show in the present work that the 6.4 keV line from this 
region is indeed most likely excited by subrelativistic ion colli- 
sions. 

The plan of the paper is as follows. In Sect. 2, we theoret- 
ically study the production of nonthermal line and continuum 
X-rays by interaction of accelerated electrons and ions with a 
neutral ambient gas. In Sect. 3, we present the XMM-Newton 



observations of the Arches cluster region and describe the data 
reduction technique we employed. In Sect. 4, we study the tem- 
poral variability of the 6.4 keV line detected from a broad region 
surrounding the star cluster. In Sect. 5, we present a detailed 
spectral analysis of the XMM-Newton data that uses the newly 
developed cosmic -ray models. The origin of the detected ther- 
mal and nonthermal radiations is discussed in Sect. 6, where we 
argue that the 6.4 keV line emission in the vicinity of the star 
cluster is produced by a large population of low-energy cosmic 
ray ions. The acceleration source of these particles is discussed 
in Sect. 7. In Sect. 8, we estimate the ionization rate induced by 
the fast ions in the ambient medium. In Sect. 9, we investigate 
the gamma-ray emission from this region. A summary is finally 
given in Sect. 10. 

2. Nonthermal X-rays from low-energy cosmic rays 

The X-ray production is calculated in the framework of a 
generic, steady-state, slab model, in which low-energy cosmic 
rays (LECRs) penetrate a cloud of neutral gas at a constant rate. 
The fast particles slow down by ionization and radiative energy 
losses in the cloud and can either stop or escape from it depend- 
ing on their path length in the ambient medium, A, which is a 
free parameter of the model. There are three other free param- 
eters that can be studied from spectral fitting of X-ray data: the 
minimum energy of the CRs entering the cloud, E^ m , the power- 
law index of the CR source energy spectrum, s, and the metal- 
licity of the X-ray emission region, Z. More details about the 
cosmic -ray interaction model are given in Appendix A. 

In Appendices B and C, we describe the atomic processes 
leading to X-ray continuum and line production as a result of 
accelerated electron and ion impacts. At this stage, we neglect 
the broad lines that can arise from atomic transitions in fast 
C and heavier ions following electron captures and excitations 
(Tatischeff et al. |19981 >. We only study the production of the nar- 
rower lines that result from K-shell vacancy production in the 
ambient atoms. We consider the Ka and K/3 lines from ambient 
C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni. We now examine the 
properties of the most important of these narrow lines in detail, 
the one at 6.4 keV from ambient Fe. 

2.1. LECR electrons 

We present in Figs. QJU characteristic properties of the X-ray 
spectrum resulting from LECR electron interactions around the 
neutral Fe Ka line. All the calculations were performed for an 
ambient medium of solar composition (i.e. Z = Z where Z Q is 
the solar metallicity). Figures Q] and [3] show the equivalent width 
(EW) and luminosity of the 6.4 keV line, whereas Figs. [2] and 
[4] show the slope of the underlying continuum emission at the 
same energy. The former two quantities depend linearly on the 
metallicity, whereas the continuum emission, which is produced 
by electron bremsstrahlung in ambient H and He, is independent 
ofZ. 

We see in Figs. QJ and [3^ that the EW of the Fe Ka line is 
generally lower than ~O.45x(Z/Z ) keV. The only exception is 
for s = 1.5 and A > 10 24 cm~ 2 . But in all cases, we expect the 
EW to be lower than 0.6 x (Z/Z ) keV. This result constitutes a 
strong constraint for a possible contribution of LECR electrons 
to the 6.4 keV line emission from the GC region, because the 
observed EW is > 1 keV in some places (see Sect. 5) and some- 
times equal to ~2 keV (see, e.g., Revnivtsev et al. 120041 ). 

The results shown in Figs. and [3^ were obtained 
without considering the additional fluorescent line emission 
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Fig. 1. Calculated a) EW and b) luminosity of the 6.4 keV Fe 
Ka line produced by LECR electrons as a function of the path 
length of the primary electrons injected in the X-ray production 
region, for five values of the electron spectral index s (Eq. ( 1A.7I )). 
The ambient medium is assumed to have a solar composition 
and the electron minimum energy E m [ n = 100 keV. In panel b, 
the luminosity calculations are normalized to a total power of 
1 erg s _1 injected by the fast primary electrons in the ambient 
medium. 




A (cm' 2 ) 

Fig. 2. Slope at 6.4 keV of the bremsstrahlung continuum emis- 
sion produced by LECR electrons, as a function of the path 
length of the primary electrons injected in the X-ray production 
region, for five values of the spectral index s. The electron min- 
imum energy is taken to be E m [ n = 100 keV. 



that can result from photoionization of ambient Fe atoms by 
bremsstrahlung X-rays >7.1 keV emitted in the cloud. This con- 
tribution can be estimated from the Monte-Carlo simulations of 
Leahy & Creighton (119931 1. who studied the X-ray spectra pro- 
duced by reprocessing of a power-law photon source surrounded 
by cold matter in spherical geometry. For the power-law photon 
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Fig. 3. Same as Fig. 1 but as a function of the electron minimum 
energy, for A = 10 25 cirT 2 . 



A = 10 25 cm" 2 




J 



10" 10 

E min (keV) 



Fig. 4. Same as Fig. 2 but as a function of the electron minimum 
energy, for A = 10 25 cm -2 . 



index a — 1 , the simulated E W of the neutral Fe Ka line can be 
satisfactorily approximated by 



EW LC 93 * 0.07 x (Z/Z ) x 



Ah 



10 23 cm 



-2 



keV 



(1) 



as long as the radial column density of the absorbing cloud, 
N£, is lower than 10 24 cm -2 . For a = 2, we have EWlc93 ~ 
0.03 x (Z/Z Q ) x (Ag/10 23 cirT 2 ) keV. Thus, we see by compar- 
ing these results with those shown in Figs. QJ and [3^ that the 
additional contribution from internal fluorescence is not strong 
for Ag < 5 x 10 23 cirT 2 . 

As shown in Figs. QJ> and[3j), the production of 6.4 keV line 
photons by LECR electron interactions is relatively inefficient: 
the radiation yield ^6.4 keV = Lx(6A ke"V)/(dW e /dt) is always 
lower than 3 x 10~ 6 (Z/Z ), which implies that a high kinetic 
power in CR electrons should generally be needed to produce 
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Fig. 5. Cross sections involved in the calculation of the Fe Ka line EW. Solid lines: cross sections (in units of barn per ambient 
H-atom) for producing the 6.4 keV line by the impact of fast electrons (left panel) and protons (right panel), assuming solar 
metallicity. Dashed lines: differential cross section (in barn per H-atom per keV) for producing 6.4 keV X-rays by bremsstrahlung 
of fast electrons (left panel) and inverse bremsstrahlung from fast protons (right panel), in a medium composed of H and He with 
H/He=0. 1 . The ratio of these two cross sections gives the EW of the 6.4 ke V line (in keV) for a mono-energetic beam of accelerated 
particles. 



an observable Ka line from neutral or low-ionized Fe atoms. 
For example, the total luminosity of the 6.4 keV line emission 
from the inner couple of hundred parsecs of our Galaxy is > 
6 x 10 34 erg s" 1 (Yusez-Zadeh et al. l20"07l l. such that dW e /dt > 
2x1 40 erg s~ 1 would be needed if this emission was entirely due 
to LECR electrons (assuming the ambient medium to be of solar 
metallicity). Such power would be comparable to that contained 
in CR protons in the entire Galaxy. 

On the other hand, Fig. QJ> shows that LECR electrons can 
produce a significant Fe Ka line (i.e. /?6.4kev ~ 10~ 6 ) in dif- 
fuse molecular clouds with iVS < 10 22 cmT 2 , especially in the 
case of strong particle diffusion for which A can be much larger 
than (see Appendix A). An observation of a 6.4 keV line 
emission from a cloud with iVjj ~ 10 21 crrT 2 would potentially 
be a promising signature of LECR electrons, since the efficiency 
of production of this line by hard X-ray irradiation of the cloud 
would be low: the ratio of the 6.4 keV line flux to the integrated 
flux in the incident X-ray continuum above 7. 1 keV (the K-edge 
of neutral Fe) is only ~ 1CT 4 for iVg = 10 21 crrT 2 (Yaqoob et al. 
[20T0T >. 

Figures |2] and |4] show the slope F of the bremsstrahlung con- 
tinuum emission, as obtained from the derivative of the differen- 
tial X-ray production rate d(dQx/dt)/dEx taken at 6.4 keV. We 
see that, for E m \ n > 100 keV, T is lower than 1.4 regardless of 
s and A. This is because bremsstrahlung X-rays < 10 keV are 
mainly produced by LECR electrons < 100 keV, and the equi- 
librium spectrum of these electrons is hard (see Fig. IB.U and 
depends only weakly on the distribution of electrons injected in 
the ambient medium at higher energies. 

Thus, after having studied the influence of the free parame- 
ters over broad ranges, we can summarize the main characteris- 
tics of the X-ray emission produced by LECR electrons as fol- 
lows. First, the continuum radiation should generally be hard, 
F < 1.4, provided that nonthermal electrons i 100 keV are 
not able to escape from their acceleration region and penetrate 
denser clouds. Secondly, the EW of the 6.4 keV Fe Ka line is 
predicted to be ~ (0.3 - 0.5) x (Z/Z Q ) keV, whatever the electron 
acceleration spectrum and transport in the ambient medium. 

The reason that the EW of the 6.4 keV line is largely inde- 
pendent of the electron energy distribution is given in Fig. [5^, 
which shows the relevant cross sections for this issue. The 
solid line is the cross section for producing 6.4 keV Fe Ka 



X-rays expressed in barn per ambient H-atom, that is ap e x 
cr^? (see Eq. ( lB.3l l). The dashed line is the differential cross 
section for producing X-rays of the same energy by electron 
bremsstrahlung. The Fe line EW produced by a given electron 
energy distribution is obtained from the ratio of the former cross 
section to the latter, convolved over that distribution. We see that 
the two cross sections have similar shapes, in particular similar 
energy thresholds, which explains why the EW of the 6.4 keV 
line depends only weakly on the electron energy. However, the 
cross section for producing the 6.4 keV line increases above 
1 MeV as a result of relativistic effects in the K-shell ioniza- 
tion process (Ouarles ll976l Kim et al. 2000b), which explains 
why the EW slightly increases with the hardness of the electron 
source spectrum (see Fig. [3^). 

Figure |5p shows the same cross sections but for proton im- 
pact. The calculation of these cross sections are presented in 
Appendix C. We see that the cross section for the line production 
has a lower energy threshold than that for the bremsstrahlung 
continuum. We thus expect that LECR protons with a relatively 
soft source spectrum can produce a higher EW of the 6.4 keV 
line than the electrons (see also Dogiel et al. 201 1 ). 



2.2. LECR ions 

Figures [6}j9] present characteristic properties of the X-ray spec- 
trum resulting from LECR ion interactions around the neutral Fe 
Ka line. As before, all the calculations were done for an ambi- 
ent medium of solar composition. The most remarkable result is 
that fast ions with a soft source spectrum can produce very high 
EW of the 6.4 keV line (Figs. [6^ and [8^). However, we see in 
Fig. [8^ that the line EW becomes almost independent of s for 
^min i 20 MeV nucleon -1 . This is because (i) the cross section 
for producing continuum X-rays at 6.4 keV and the one for the 
neutral Fe Ka line have similar shapes above 20 MeV nucleon 1 
(see Fig. |5p) and (ii) the CR equilibrium spectrum below 
only weakly depends on s and A. 

As shown in Fig.|6j5, for E m i n = 1 MeV nucleon 1 , the ra- 
diation yield Re.4k&v can reach ~10~ 6 only for relatively hard 
source spectra with s < 2 and for A i 10 25 crrT 2 , which should 
generally mean strong particle diffusion in the X-ray produc- 
tion region. For such a high CR path length, the X-rays are 
mainly produced in thick- target interactions. An efficiency of 
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Fig. 6. Calculated a) EW and b) luminosity of the 6.4 keV 
Fe Ka line produced by LECR ions as a function of the path 
length of the CRs injected in the X-ray production region, for 
five values of the CR source spectral index s (Eq. ( lA.7t ). The 
ambient medium is assumed to have a solar composition, and 
the minimum energy of the CRs that penetrate this medium is 
Emm = 1 MeV nucleon -1 . In panel b, the luminosity calcula- 
tions are normalized to a total power of 1 erg s -1 injected by the 
fast primary protons in the ambient medium. 
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Fig. 7. Slope at 6.4 keV of the continuum emission produced 
by LECR ions as a function of the path length of the fast ions 
in the X-ray production region, for five values of the spectral 
index s. The minimum energy of injection is taken to be Emm = 
1 MeV nucleon 1 . 



~1Q 6 in the production of the 6.4 keV line can also be achieved 
with a softer CR source spectrum, if E m [ n is in the range 10- 
100 MeV nucleon -1 (Fig. [8}?). This is because then most of the 
CRs are injected into the X-ray production region at energies 
where the cross section for producing Fe Ka X-rays is highest 
(see Fig. |5]3). But in any case, we find that to get 7?6.4 keV ~ 10 
it requires A > 10 24 cm -2 (for solar metallicity). It is another 
difference from the production of nonthermal X-rays by LECR 
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Fig. 8. Same as Fig. 6 but as a function of the minimum energy 
of injection, for A = 10 25 cm 2 . 
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Fig. 9. Same as Fig. 7 but as a function of the minimum energy 
of injection, for A = 10 25 crrT 2 . 



electrons, for which ^64keV can reach ~10 6 for A as low as 
10 22 cm" 2 (Fig.HJ)). 

Figures|7]and|9]show that the characteristic power-law slope 
of the continuum emission around 6.4 keV can vary from ~1 
to ~6 for s in the range 1.5-5. However, for s < 2, which is 
expected for strong shock acceleration of nonrelativistic parti- 
cles, and A > 10 24 cirT 2 , which can result from strong particle 
diffusion in the cloud, we expect F between 1.3 and 2. For these 
conditions, the EW of the neutral Fe Ka line is predicted to be in 
the naiTow range (0.6 - 0.8) x (Z/Z Q ) keV (Figs.|6ji and[8ji). This 
result could account for EW i 1 keV as observed from regions 
near the GC, provided that the diffuse gas there has a super-solar 
metallicity (i.e. Z > Z Q ). 
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Table 1. Summary of the XMM-Newton/EPIC observations available for the Arches cluster. 



Date 


Observation 
ID 


Good time exposure 
(ks) 


MOS 
Noisy CCD 






covera 


ge (%) 






Cloud 


Cluster 


Ml 


M2 


pn 


Ml 


M2 


pn 


Ml 


M2 


pn 


2000-09-19 


0112970401 


21 


21 


16 


- 


100 


100 


86 


100 


100 


100 


2000-09-21 


0112970501 


7 


9 




3 


- 


96 


99 


87 


99 


100 


100 


2001-09-04 


0112972101 


20 


20 


17 


- 


100 


100 


40 


100 


100 


21 


2002-02-26 


0111350101 


43 


42 


36 


- 


100 


100 


56 


100 


100 


94 


2002-10-02 


0111350301 


6 


6 




3 


- 


99 


100 


47 


100 


100 


82 


2004-03-28 


0202670501 


10 


14 




2 


- 


100 


100 


100 


100 


100 


100 


2004-03-30 


0202670601 


27 


27 


17 


- 


100 


100 


99 


100 


100 


100 


2004-08-31 


0202670701 


64 


74 


40 


M2-5 


100 


100 


100 


100 


100 


100 


2004-09-02 


0202670801 


83 


89 


50 


M2-5 


100 


100 


100 


100 


100 


100 


2006-02-27 


0302882601 


2 


2 




1 


- 





100 


55 





100 


91 


2006-09-08 


0302884001 


6 


6 




4 


Ml-4 


100 


99 


33 


100 


100 


7 


2007-02-27 


0506291201 


18 


20 




a 


- 


100 


87 


a 


100 


96 


a 


2007-03-30 


0402430701 


23 


24 


16 


M2-5 





100 


100 





100 


100 


2007-04-01 


0402430301 


54 


58 


29 


Ml-4, M2-5 





100 


99 





100 


100 


2007-04-03 


0402430401 


39 


40 


23 


Ml-4 





100 


100 





100 


100 


2007-09-06 


0504940201 


8 


9 




5 


M2-5 


98 


99 


34 


98 


100 


23 


2008-03-04 


0511000301 


4 


4 




2 







100 


64 





100 


100 


2008-03-23 


0505670101 


68 


73 


51 


Ml-4, M2-5 





100 


b 





100 


b 


2008-09-23 


0511000401 


4 


4 




4 


Ml-4, M2-5 


99 


99 


41 


98 


100 


41 


2009-04-01 


0554750401 


32 


32 


24 







100 


98 





100 


100 


2009-04-03 


0554750501 


40 


41 


31 


Ml-4 





100 


99 





100 


100 


2009-04-05 


0554750601 


36 


37 


25 


Ml-4 





100 


99 





100 


100 


Total 

exposure (ks) 


Imaging 


615 


652 


399 




Spectroscopy 


317 


652 


Cloud 


Cluster 


276 


315 



Notes. For each observation and instrument, we report the total good-time interval exposure after flare screening. For a number of observations, 
the Arches cluster lies at the position of the CCD 6 of MOS 1, which is out of service. We indicate the CCDs identified as noisy for MOS (data 
not used for the analysis) and the spatial coverage of the pn (incomplete due to dead columns) for our two spectral extraction regions (see Fig.llOb. 
Only observations with a spatial coverage of the region greater than 85% were used for the spectral analysis. 

{a) The pn camera was in timing mode for this observation. {b) The pn data of this observation were not taken into account for the spectrum 
extraction, because the pn spatial coverage of the background region CFig.l 10b was to low. 



3. XMM-Newton observations and data reduction 

3.1. Data reduction and particle background substraction 

For our analysis, we have considered all public XMM- 
Newton EPIC observations encompassing the Arches cluster 
(R.A.= 17 h 45 m 50 s , Dec=-28°49'20"). The criteria were to have 
more than 1 .5 ks observation time available for each camera, to 
be in full frame or extended full frame mode, and to use the 
medium filter. The data was reduced using the SAS software 
package, version 10.0. Calibrated-event files were produced us- 
ing the tasks EMCHAIN for the MOS cameras and EPPROC 
for the pn camera. We excluded from the analysis the period 
contaminated by soft proton flares, by using an automatic 3o~- 
clipping method (Pratt & Arnaud 2 003). Table[TJprovides the list 
of the 22 selected observations and their respective observing 
time per instrument after flare rejection. 

We searched for any anomalous state of MOS CCD chips 
(Kuntz & Snowden 2008 ) by performing a systematic inspection 
of the images and spectra of each chip in the 0.3-1 keV energy 
band. We identified 14 occurrences of a noisy chip in the list of 
observations (see Table Q]). The chips affected in our observa- 
tions by a high-level, low-energy background state are CCD 4 of 
MOS 1 and CCD 5 of MOS 2. We excluded data of those chips 
from our analysis when they were noisy. 

For MOS cameras, we selected events with PATTERN < 12. 
Only events with PATTERN < 4 were kept for the pn instrument. 



Depending on the nature of the analysis, we defined two kinds 
of quality-flag selection: 

- For imaging, to select events with good angular reconstruc- 
tion, we used the flags XMMEAJ3M and XMMEA.EP for the 
MOS and pn cameras, respectively. 

- For the spectrum analysis, we chose events with 
FLAG=XMMEA_SM for MOS cameras (good energy 
reconstruction) and with FLAG=0 for the pn camera. 

The particle background was derived from filter-wheel 
closed (FWC) observations that were compiled until revolution 
about 1600. To be consistent, we applied exactly the same event 
selection criteria to both data and FWC files. We checked that 
even if the particle flux has increased significantly between 2000 
and 2009, its spectrum and spatial repartition have not signifi- 
cantly changed during that period. We used the count rates be- 
tween 10 and 12 keV for MOS and between 12 and 14 keV 
for pn to normalize the FWC background level to that of our 
observations. Regions with bright sources in the observations 
have been excluded to calculate the normalization factor. Finally, 
the EVIGWEIGHT task was used to correct vignetting effects 
(IPratt et al. 20071 1. 

3.2. Maps generation 

For each observation and instrument, we produced count im- 
ages (EVSELECT task) in two energy bands (6.3-6.48 keV and 
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Fig. 10. XMM-Newton/EPIC continuum-subtracted Fe Ka emission line maps of the Arches cluster region at 6.4 keV (left panel) and 
6.7 keV (right panel). The images have been adaptively smoothed at a signal-to-noise of 20. The magenta circle indicates the region 
("Cluster") used to characterize the Arches cluster X-ray emission, which shows strong Fe Ka emission at 6.7 keV. The region 
inside the white ellipse but outside the magenta circle indicates the region ("Cloud") used for spectral extraction to characterize the 
bright 6.4 keV regions surrounding the Arches cluster. The region inside the yellow circle but outside the two dashed ellipses shows 
the local background used for the spectral analysis. The axes of the maps (in green) indicate Galactic coordinates in degrees. North 
is up and east to the left. 



6.564-6.753 keV), which are dominated by the Fe Ka lines at 
6.4 keV from neutral to low-ionized atoms and at 6.7 keV from 
a hot thermally-ionized plasma. For each energy band, observa- 
tion, and instrument, the normalized particle background image 
derived from FWC observations was subtracted from the count 
image. The particle background events were rotated beforehand 
so as to match the orientation of each instrument for each ob- 
servation. For each energy band, observation, and instrument, an 
exposure map was generated (EEXMAP task) taking the different 
efficiencies of each instrument into account. 

To produce line images of the Fe Ka emission at 6.4 and 
6.7 keV, the continuum under the line needs to be subtracted. 
For that, we produced a map in the 4. 17-5.86 keV energy band, 
which is dominated by the continuum emission, with the same 
procedure as before. To determine the spectral shape of the con- 
tinuum in the 4-7 keV band, the spectrum of the cloud region 
(representative region defined in Fig. [10] and Table |2]l was fitted 
by a power law and Gaussian functions to account for the main 
emission lines. The continuum map was then renormalized to 
the power-law flux in the considered energy band and subtracted 
from the corresponding energy band image. 

The background-subtracted and continuum-subtracted count 
images and the exposure maps were then merged using the 
EMOSAIC task to produce images of the entire observation set. 
The resulting count images were adaptively smoothed using the 
ASMOOTH task with a signal-to-noise ratio of 20. The count 
rate images were obtained by dividing the smoothed count im- 
ages by the smoothed associated exposure map. The template of 
the smoothing of the count image was applied to the associated 
merged exposure map. 

Figure [10] shows the resulting Fe Ka line maps at 6.4 and 
6.7 keV of the Arches cluster region. The star cluster exhibits 
strong Fe Ka emission at 6.7 keV. Bright Fe Ka 6.4 keV struc- 
tures are observed around the Arches cluster. 



Table 2. Definition of the spectral extraction regions. 



Region 


RA (J2000) 


Dec (J2000) 


Shape 


Parameters 


Cluster 


17 h 45 ra 50.3 s 


-28°49'19" 


circle 


15" 


Cloud 
excl. 


17 h 45 m 51.0 s 
17 h 45 ra 50.3 s 


-28°49'16" 
-28°49'19" 


ellipse 
circle 


25", 59", 155° 
15" 


Bkg 
excl. 
excl. 


17 h 45 ra 51.0 s 
17 h 45 m 47.2 s 
17 h 45 ra 50.4 s 


-28°49'25" 
-28°50'42" 
-28°49'03" 


circle 

ellipse 

ellipse 


148" 

37", 78", 130° 
55", 100", 160° 



Notes. This Table provides the center position, circle radius and for 
ellipses, minor, and major axes, as well as angle (counter clockwise 
from straight up). "Bkg" means background and "excl." indicates the 
zones of exclusion. 



3.3. Spectrum extraction 

To characterize the properties of the emission of the Arches 
cluster and its surroundings, we defined two regions, which are 
shown in Fig.[l0] The region called "Cluster" corresponds to the 
Arches cluster and exhibits strong Fe Ka emission at 6.7 keV. 
The region called "Cloud" corresponds to the bright Fe Ka 6.4 
keV emission structures surrounding the Arches cluster. Table [2] 
provides the coordinates of these regions. 

In 9 of the 22 relevant observations, the Arches region lies 
on the out-of-service CCD6 of MOS 1, reducing the available 
observation time by more than a factor 2 compared to MOS 2 
(see Table [TJ. Regarding the pn camera, the cloud and cluster 
regions suffer from the presence of dead columns in a number 
of observations. We thus estimated the spatial coverage of each 
region for each pn observation (see TableQ}. After several tests, 
we chose for the spectral analysis to keep only observations with 
a pn spatial coverage greater than 85%. For all these selected ob- 
servations, the MOS spatial coverage of each region was greater 
than 85% (see Table Q}. 
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Fig. 11. Lightcurve of the 6.4 keV Fe Ka line flux arising from 
a large region around the Arches cluster (the region labeled 
"Cloud" in Fig. ITOb . The red horizontal line shows the best fit 
with a constant flux. 



Table Q] summarizes the final total available exposure time 
by instrument for spectral analysis. With the 22 observations, 
we obtained 317 ks for MOS 1 and 652 ks for MOS 2. For the 
pn spectral analysis, we obtained 276 and 315 ks on the cloud 
and cluster regions, respectively. 

For each region, the particle background spectrum was esti- 
mated from the FWC observations in the same detector region. 
The astrophysical background around the Arches cluster shows 
spatial structures and, notably, an increase towards the Galactic 
plane. After several tests, we concluded that the most represen- 
tative local background for the Arches region was that of the 
region encircling the Arches, but avoiding the zones emitting 
at 6.4 keV. The background region is defined in Table [2] and 
shown in Fig. \W\ (yellow circle with the exclusion of the two 
dashed-line ellipses). To subtract the particle and astrophysical 
background from the spectra, we used the method of double sub- 
traction described in |Arnaud~e t al. (2002), The ancillary and re- 
distribution matrix function response files were generated with 
the SAS ARFGEN and RMFGEN tasks, respectively. 

The spectra from individual observations of the same region 
were then merged for each instrument and rebinned to achieve a 
signal-to-noise per bin of 3 cr. 

4. Variability of the 6.4 keV line 

The temporal variability of the 6.4 keV line flux is a key di- 
agnostic for deciphering the origin of the line (see, for exam- 
ple, Ponti et al. 1201 Oi l. To study this aspect, we combined spec- 
tra extracted from the cloud region to obtain a sampling of the 
emission at seven epochs: September 2000 (2 observations), 
February 2002 (1 observation), March 2004 (2 observations), 
August/September 2004 (2 observations), March/ April 2007 (3 
observations), March 2008 (2 observations), and April 2009 (3 
observations). This sampling is similar to the one recently used 
by Capelli et al. (1201 lbk except that, for an unknown reason, 
these authors did not included the September 2000 epoch in their 
analysis. 

To measure the intensity of the neutral or low-ionization Fe 
Ka line at each epoch, we modeled the X-ray emission from the 
cloud region as the sum of an optically thin, ionization equilib- 
rium plasma (APEC, Smith et al. 2001 ), a power-law continuum, 
and a Gaussian line at ~6.4 keV. These three components were 



subject to a line-of-sight photoelectric absorption so as to ac- 
count for the high column density of the foreground material. 
We used the X-ray spectral-fitting program XSPEC3 to fit this 
model simultaneously to EPIC MOS and pn spectra between 1 .5 
and 10 keV. More details on the fitting procedure will be given in 
the next section. All the fits were satisfactory and gave reduced 
X 1 ~ 1. 

The photon fluxes in the 6.4 keV line thus determined are 
shown in Fig. QT| The best fit of a constant flux to these data 
is satisfactory, giving a x 2 of 3.3 for six degrees of freedom 
(dof). The significance of a variation in the line flux is then 
only of 0.3oo The best-fit mean flux is F 6.4 keV = (8.8 + 0.5) x 
10~ 6 ph cirT 2 s" 1 . 

The fact that the intensity of the 6.4 keV line emitted from 
the vicinity of the Arches cluster is consistent with being con- 
stant is in good agreement with the previous work of Capelli et 
al. (1201 lbl . We note, however, that the line fluxes obtained in 
the present work are systematically lower by ~20%. This is at- 
tributable to a difference in the background modeling: whereas 
we used a broad region as close as possible to the Arches clus- 
ter to subtract the astrophysical background prior to the spectral 
fitting (Fig. [Toll. Capelli et al. included the background as a com- 
ponent within the fitting. 



5. X-ray spectral analysis of time-averaged spectra 

We again used XSPEC to fit various models to time-averaged 
spectra extracted from the two source regions shown in Fig.fTOl 
The fits were performed simultaneously on the stacked MOS1, 
MOS2, and pn spectra, but we allowed for a variable cross- 
normalization factor between the MOS and pn data. Independent 
of the fitting model, we found very good agreement between the 
MOS and pn cameras, to better than 1 % for the data extracted 
from the cluster region and to 4-5% for the data extracted from 
the cloud region. These factors are consistent with the resid- 
ual uncertainty in the flux cross-calibration of the EPIC cameras 
(Mateos et al. 120091 . 



5. 1 . X-ray emission from the star cluster 

We first modeled the emission of the cluster region as the sum 
of an APEC plasma component and a nonthermal component 
represented by a power-law continuum and a Gaussian line at 
~6.4 keV. The centroid energy of the Gaussian line was allowed 
to vary, but the line width was fixed at 10 eV. All the emis- 
sion components were subject to a line-of-sight photoelectric ab- 
sorption (WABS model in XSPEC). The best-fit results obtained 
with this model (called model 1 in the following) are reported in 
Tableland the corresponding spectra shown in Fig.[T2"k In this 
table and in the following discussion, all the quoted errors are at 
the 90% confidence level. 

This fitting procedure did not allow us to reliably constrain 
the metallicity of the X-ray emitting plasma. Indeed, the best 
fit was obtained for a super-solar metallicity Z > 5 Z , which 
is not supported by other observations. Such an issue has al- 
ready been faced in previous analyses of the X-ray emission 



1 http : //heasarc . nasa . gov/xanadu/ xspec/ 

2 Noteworthy is that these results were obtained without 
taking the systematic error in the effective area of the EPIC 
camera into account. This error is estimated to be 7% for 
on-axis sources and to increase with off-axis angle (see the 
XMM-Newton Calibration Technical Note CAL-TN-0018.pdf at 



http : / /xitim . vilspa . esa . es/external/xmm_sw_cal/calib/ documentat 
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Table 3. Spectral analysis of the X-ray emission from the Arches star cluster and associated cloud region with standard XSPEC 
models. 





(Unit) 




star cluster 




cloud region 






model 1 


model 2 


model 3 


model 1 




y ru 1 1 cm ) 






Q Q+O.St 




ft. i yz.) 


VJ^c v ) 




27+0 04 


92 +0 14 




z/z e 




- 


1.7 (fixed) 


1.7 (fixed) 


- 


Lt(2) 


(see notes below) 




1100+ 1500 


13 +7 




Nu(l) 


(10 22 Ff crrr 2 ) 


9.5+0.3 


12 +ub 


12 8 +l j 


11 3 +w 


kT(l) 

V A ) 

z/z e 


(keV) 


i 70+U.O6 
• -0.05 

1.7 (fixed) 


1,u -0.05 

1.7 (fixed) 


1 78 +U15 
1.7 (fixed) 


9 9+1.0 

^ -0.5 

1.7 (fixed) 


4r(l) 


(see notes below) 


20.4±1.8 


30±4 


23+4 


5 2 +4 - 9 


E(,A keV 


(keV) 


6.41+0.02 


6.40+0.02 


6.40±0.02 


6.409+0.005 


F 6.4 kcV 


(1(T 6 phcrrT 2 s" 1 ) 


1.2±0.3 


1.3+0.3 


1.3±0.3 


o 9+0.5 


r 

V 


(lfr 5 cm" 2 s-» keV"') 


0.7±0.4 
1 3 +1 - 4 


4 +u 5 
62 +1 15 


1 3 +2S 

^ -0.9 


1 -"-Q.2 


EW 6 . 4 kcV 


(keV) 


0.4±0.1 


0.4±0.1 


0.4±0.1 


1.2±0.2 


^/dof 




1222/978 


1152/976 


1129/975 


560/491 



Notes. Model 1: WABSx(APEC + Gaussian + powerlaw); model 2: WABSx(APEC + APEC + Gaussian + powerlaw); model 3: WABSxAPEC 
+ WABSx(APEC + Gaussian + powerlaw). A^: absorption column density. kT, Z/Z , and I/, T : temperature, metallicity relative to solar, and 
normalization of the APEC thermal plasma (I kT is in unit of 10~ 18 J n^n-^dVKAixD 1 ), where n e and n H are the electron and proton number densities 
(cm -3 ) and D the distance to the source in cm). E^a keV and F^a ^v'- centroid energy and flux of the neutral or low-ionization Fe Ka line. T and 
/p i : index and normalization at 1 keV of the power-law component. EW 6 4 ke v: EW of the 6.4 keV line with respect to the power-law continuum. 
^ 2 /dof: x 1 P er degree of freedom. 




Energy (keV) Energy (keV) 

Fig. 12. X-ray spectra of the Arches cluster as measured in the MOS 1 (black), MOS2 (red), and pn (green) cameras aboard XMM- 
Newton, compared to a) a model with only one thermal plasma component (model 1 in Table [3} and b) a model with two thermal 
plasma components (model 3 in Table[3]l. The lower panels show the associated residuals in terms of standard deviations. The second 
plasma of temperature kT = 0.9 keV accounts for a significant emission in the He-like Si and S Ka lines at 1.86 and 2.46 keV, 
respectively. 



from the Arches cluster region. Thus, Tsujimoto et al. (120071 1 
fixed the plasma metallicity to be solar in their analysis of Suzaku 
data, whereas Capelli et al. ( 1201 lal 1201 lbl ) adopted Z = 2 Z 
in their analysis of XMM-Newton data. With the Chandra X-ray 
Observatory, Wang et al. (2006) were able to resolve three bright 
point-like X-ray sources in the core of the Arches cluster - most 
likely colliding stellar wind binaries - and study them individ- 
ually. These sources were all modeled by an optically thin ther- 
mal plasma with a temperature of ~1.8— 2.5 keV and a metal- 
licity Z/Z = 1-8+Q 2- In our analysis, we fixed the metallicity 
of the thermal plasma in the cluster region to be 1 .7 Z Q , which 
is the best-fit value that we were able to obtain for the cloud 



region using the LECR ion model developed in this paper to ac- 
count for the nonthermal emission (see Sect. 5.2 below). The 
adopted metallicity is also consistent with the results of Wang et 
al. d2006t . 

Model 1 gives a good fit to the data from the cluster region 
above ~3 keV. In particular, the detection of the neutral or low- 
ionization Fe Ka line is significant (see Table 0). But the fit is 
poorer below 3 keV, because the data shows clear excesses of 
counts above the model at ~1.85 and ~2.45 keV (see Fig.[T2h). 
These features most likely correspond to the Ka lines from He- 
like Si and S, respectively. We checked that this excess emission 
is not due to an incomplete background subtraction by produc- 
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Fig. 13. Same as Fig. [TO] but for the He-like Si Ka line at 
1.86 keV. 



ing a Si Ka line image in the energy band 1.76-1.94 keV. To 
estimate the contribution of the continuum under the Si line, we 
first produced a count map in the adjacent energy band 2.05- 
2.15 keV and then normalized it to the expected number of con- 
tinuum photons in the former energy range. The normalization 
factor was obtained from a fit to the EPIC spectra of the cluster 
region by model 1 plus two Gaussian functions to account for 
the Si and S Ka lines. The resulting map in the Si line shows 
significant excess emission at the position of the Arches cluster 
(FigEJ. 

To account for the presence of the He-like Si and S Ka lines 
in the X-ray spectrum of the cluster, we included in the fit- 
ting model a second APEC component subject to the same 
photoelectric absorption as the other components (model 2). 
The quality of the fit significantly improves with this additional 
thermal plasma component (y 2 =1152 for 976 degrees of free- 
dom), whose best-fit temperature is kT - 0.27 + 0.04 keV 
(Table [3}. However, the absorption-corrected intrinsic luminos- 
ity of this plasma is found to be quite high in the soft X-ray 
range: Lj nt (0.4 - 1 keV) « 2.3 x 10 36 erg s _1 , assuming the dis- 
tance to the GC to be D = 8 kpc (Ghez et al. 120081 . 

In a third model, we let the X-ray emission from the plasma 
of lower temperature be absorbed by a different column den- 
sity than the one absorbing the X-rays emitted from the other 
components. It allows us to further improve the fit to the data 
(X 2 -\\29 for 975 degrees of freedom, Table|3j see also Fig.fT2b 
for a comparison of this model to the data). We then found 
kT = 0.92+°;^ keV and L int (0.4 - 1 keV) * 2.0 x 10 34 erg s" 1 for 
the lower temperature plasma. The origin of this thermal com- 
ponent, which was not detected in previous X-ray observations 
of the Arches cluster, is discussed in Sect. 6.1 below. 

As can be seen from Table[3] the addition of a second APEC 
component in the fitting model of the star cluster emission in- 
creases the absorbing column density Nn(Y) significantly. It also 
has some impact on the temperature of the hotter plasma and on 
the index of the power-law component, but not on the properties 
of the 6.4 keV line. 

5.2. X-ray emission from the cloud region 

We used model 1 to characterize the X-ray emission from the 
cloud region, except that we also included a Gaussian line at 
7.05 keV (fixed centroid energy) to account for the neutral or 
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Table 4. Spectral analysis of the X-ray emission from the cloud 
region with LECR electron and ion models. 





(Unit) 


T pr^R pIpcItotiq 


T FCR ions 


A! 


(1U xi cm ) 


11.9_ 14 


11 0+1,4 

12.2_, 


kT 

z/z a 

hr 


(keV) 
(see notes below) 


1 q+0.6 

> 3.1 


2 +u ? 

^•"-0.3 

1.7±0.2 
7.0^-° 


A 

s 

F 

^LECR 


(H-atoms cirT 2 ) 

(keV) or (keV/n) 
(10~ 8 erg cirr 2 S" 1 ) 


5 x 10 24 (fixed) 
> 2.5 
< 41 
5 +74 


5 x 10 24 (fixed) 

j n+0.5 
-0.6 

10 4 (fixed) 

5 6 +0 - 7 


*7dcf 




558/492 


558/493 



Notes. XSPEC model: WABSx(APEC + LECR//), where p stands for 
electrons or ions. 7V H , kT, Z/Z e , and I kT : as in Table[3] A, s, E min , and 
A'lecr: LECR path length, source spectrum index, minimum energy, 
and model normalization. By definition dW/dt = 4nD 2 N LEC R is the 
power injected in the interaction region by primary CR electrons or pro- 
tons of energies between E min and £ max = 1 GeV (D is the distance to 
the source). 



low-ionization Fe K/3 line. The Fe K/?/Ka flux ratio was imposed 
to be equal to 0.13 (Kaastra & Mewe 119931 . We checked that 
including a second thermal plasma component (model 2 or 3) 
is not required for this region, as it does not improve the qual- 
ity of the fit. As before, we fixed the metallicity of the emitting 
plasma to be 1.7 Z Q . The best-fit temperature, kT = 2.2^" keV, 
is marginally higher than the one of the high-temperature plasma 
emanating from the cluster region. 

We now compare the characteristics of the prominent non- 
thermal emission of the cloud region with the model predictions 
discussed in Sect. 2. In the LECR electron model, the measured 
value r = 1.6jj would only be expected for low values of 
the CR minimum energy E m i n < 100 keV and for relatively soft 
source spectra with s i 2.5 (see Fig. HJ. But for these CR spec- 
trum parameters the neutral Fe Ka line is predicted to be rela- 
tively weak, EW 6 . 4 ke v < 0.4 x (Z/Z Q ) keV (see Fig. 0. Thus, 
it would require an ambient Fe abundance i 3 times the solar 
value to account for the measured EW of 1.2 + 0.2 keV. The 
measured properties of the nonthermal component emitted from 
the cloud region thus appear to be hardly compatible with the 
predictions of the LECR electron model. 

On the other hand, the measured values of T and EW6.4 kev 
for the cloud region seem to be compatible with the LECR ion 
model. The measured power-law slope can be produced in this 
model with any spectral index s ~ 1 .5-2, provided that the CR 
path length A > 10 24 cm -2 (Figs. Q and |9). We then expect 
EW 6 . 4 kev ~ (0.6-1) x (Z/Z ) keV (Figs. [6] and ©, which would 
be in good agreement with the measured EW for an ambient 
metallicity Z < 2Z . 

To further study the origin of the prominent nonthermal 
emission of the cloud region, we created LECR electron and ion 
models that can be used in the XSPEC software. For this pur- 
pose, a total of 70875 spectra were calculated for each model 
by varying the four free parameters of the models in reasonable 
ranges. The calculated spectra were then gathered in two FITS 
files that can be included as external models in XSPECfl We 
then fitted the stacked spectra of the cloud region by the XSPEC 
model WABSx(APEC + LECR/?), where p stands for electrons 
or ions. The best-fit results obtained with both models are given 
in Table |4] 

3 These models are available upon request to the authors. 
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Energy (keV) Energy (keV) 



Fig. 14. a) X-ray spectra of the cloud region as measured in the XMM-Newton cameras and the best-fit spectral model assuming 
that the emission comes from a combination of a collisionally ionization equilibrium plasma (APEC model) and a nonthermal 
component produced by interactions of LECR ions with the cloud constituents (see Table[4]l; b) model components. 



In this spectral fitting, we allowed for a variable metallicity 
of the nonthermal X-ray production region (i.e. the parameter 
Z of the LECRp models), but we imposed this parameter to be 
equal to the metallicity of the thermal plasma. Since both fits did 
not usefully constrain the path length of the LECRs in the in- 
teraction region, we fixed A = 5 x 10 24 cirT 2 for both models, 
which, as discussed in Appendix A, is a typical value for nonrel- 
ativistic protons propagating in massive molecular clouds of the 
GC environment (see Eq. (IA.6I 1). As anticipated, the LECR elec- 
tron model cannot satisfactorily account for the data, because the 
best fit is obtained for too high a metallicity (Z > 3.1 Z ; limit 
at the 90% confidence level) and a low CR minimum energy 
(E m m < 41 keV). This conclusion is independent of the adopted 
value of A. 

On the other hand, the data can be characterized well by a 
thin plasma component plus an LECR ion model. In particu- 
lar, the best-fit metal abundance Z/Z — 1.7 + 0.2 is in good 
agreement with previous works (Wang et al. [2006). The best- 
fit CR spectral index is s = 1.9*2*. For such a relatively hard 
CR source spectrum, one can see from Figs. [8^ and [9] that 
the nonthermal X-ray emission produced by LECR ions only 
weakly depends on the CR minimum energy Emit,. Accordingly, 
the fit did not constrain this parameter, which was finally fixed 
at E mm = 10 MeV nucleon -1 . As discussed in Appendix A, 
the process of CR penetration into molecular clouds is not un- 
derstood well, such that Emm is loosely constrained from the- 
ory. This parameter has an effect, however, on the power in- 
jected by the primary CRs in the X-ray production region (see 
Fig. |8p). Thus, for E m \ n = 1 MeV nucleon 1 (resp. E m m = 
100 MeV nucleon -1 ), the best-fit normalization of the LECR 
ion model is A^lecr = (7.4+jg) x 10~ 8 erg cirT 2 s -1 (resp. 
(3.1^q5)x 10 -8 erg cirT 2 s ). The corresponding power injected 
by LECR protons in the cloud region (dW/dt = 47rD 2 A^LECR) lies 
in the range (0.2-l)xl0 39 erg s _1 (with D = 8 kpc). 

The best-fit model obtained with the LECR ion component 
is compared to the data of the cloud region in Fig. [14k and the 
corresponding theroretical spectrum is shown in Fig. [T4b . The 
latter figure exhibits numerous lines arising from both neutral 
and highly-ionized species, which could be revealed by a future 
instrument having an excellent sensitivity and energy resolution. 



6. Origin of the detected radiations 

We have identified three distinct components in the X-ray spec- 
tra extracted from the cluster region: an optically thin thermal 
plasma with a temperature fcr~1.6-1.8 keV, another plasma of 
lower temperature (kT ~0.3 keV in model 2 or ~0.9 keV in model 
3), and a relatively weak nonthermal component characterized 
by a hard continuum emission and a line at 6.4 keV from neutral 
to low-ionized Fe atoms (EW6.4 keV = 0.4 + 0.1 keV). The X- 
ray radiation arising from the cloud region is also composed of a 
mix of a thermal and a nonthermal component, but the 6.4 keV 
Fe Ka line is much more intense from there, with a measured 
EWof 1.2 + 0.2keV. 

6.1. Origin of the thermal X-ray emissions 

The thermal component of temperature W~1.6-1.8 keV de- 
tected from the star cluster most likely arises from several 
colliding stellar wind binaries plus the diffuse hot plasma of 
the so-called cluster wind. Wang et al. (2006) find with the 
Chandra telescope three point-like sources of thermal emission 
with kT~l. 8-2.5 keV embedded in a spatially extended emis- 
sion of similar temperature. Capelli et al. (I201 lal l have recently 
found with XMM-Newton that the bulk of the X-ray emission 
from the Arches cluster can be attributed to an optically thin ther- 
mal plasma with a temperature kT~\.l keV. The diffuse thermal 
emission from the cluster is thought to be produced by the ther- 
malization of massive star winds that merge and expand together. 
The expected temperature of such a cluster wind is consistent 
with the temperature of the hot thermal component identified in 
this and previous works (see Capelli et al. 201 la and references 
therein). 

The plasma with kT~\.6— 1.8 keV is at the origin of the He- 
like Fe Ka line at 6.7 keV. The corresponding map generated 
in the present work (Fig. [10] right panel) is in good agreement 
with the Chandra observations. Wang et al. (2006) suggest that 
the observed elongation of this emission in the east-west direc- 
tion reflects an ongoing collision of the Arches cluster with a 
local molecular cloud traced by the CS emission. As discussed 
by Wang et al. (2006), this collision may help in explaining the 
spatial confinement of this hot plasma. 
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The second thermal component of temperature kT~0.3 keV 
(model 2) or kT~0.9 keV (model 3) was not detected in previ- 
ous X-ray observations of the Arches cluster. An optically thin 
thermal plasma of temperature IcT~0.8 keV was reported by 
Yusef-Zadeh et al. (2002b) from Chandra observations, but not 
by subsequent X-ray observers (Wang et al. 2006 , Tsujimoto et 
al. 120071 Capelli et al. 1201 la fr. However, the high-quality spec- 
tral data obtained in the present work reveal that a single APEC 
thermal plasma model cannot account simultaneously for the ob- 
served lines at ~ 1 .85, ~2.45, and 6.7 keV, which arise from He- 
like Si, S, and Fe atoms, respectively. The map at ~ 1 .85 keV 
clearly shows that the star cluster significantly emits at this en- 
ergy (Fig. [13]). The required additional plasma component is sub- 
ject to a high interstellar absorption: Nh ~ 1.2 x 10 23 in model 
2 and 8.3 x 10 22 H cirT 2 in model 3 (Table [3]). It shows that the 
emitting plasma is located in the Galactic center region and not 
in the foreground. 

The temperature of the second thermal component suggests 
that this emission could be due to a collection of individual mas- 
sive stars in the cluster. Single hot stars with spectral types O and 
early B are known to emit significant amounts of thermal X-rays 
with a temperature kT in the range 0.1-1 keV and a typical lumi- 
nosity in soft X-rays L x (0A- 1 keV) ~ 1.5 x 10~ 7 Lboi (Antokhin 
et al. 120081 Giidel & Naze |2"0"09j l. Here, L boX is the bolometric 
luminosity of the star. The total bolometric luminosity of the 
Arches cluster is ~ 10 7 ° L Q and most of it is contributed by 
early B- and O-type stars, some of which have already evolved 
to the earliest Wolf-Rayet phases (Figer et al. 2002). Then, a to- 
tal soft X-ray luminosity L X (0A - 1 keV) ~ 3.6 x 10 34 erg s" 1 
can be expected from the ensemble of hot massive stars of the 
cluster. This estimate is much lower than the unabsorbed in- 
trinsic luminosity of the ~0.3 keV plasma found in model 2, 
£int(0.4 - 1 keV) * 2.3 x 10 36 erg s" 1 . But it is roughly con- 
sistent with the absorption-corrected luminosity of the ~0.9 keV 
plasma found in model 3: L int (0.4- 1 keV) * 2.0x 10 34 erg s It 
is not clear, however, why the latter component is less absorbed 
than the high-temperature plasma emitted from the Arches clus- 
ter (Table©. 



6.2. Origin of the 6.4 keV line emission 
6.2.1. The cloud region 

Several molecular clouds of the GC region emit the 6.4 keV line, 
most notably Sgr Bl, Sgr B2, Sgr C, and clouds located be- 
tween Sgr A* and the Radio Arc (see Yusef-Zadeh et al. 2007 ). 
Detections of time variability of the 6.4 keV line from Sgr B2 
(Inui et al. I2009I I, as well as from molecular clouds within 15' 
to the east of Sgr A* (Muno et al. [20071 Ponti et al. EoTUl i, are 
best explained by the assumption that the Fe Ka line emission 
from these regions is a fluorescence radiation produced by the 
reprocessing of a past X-ray flare from the supermassive black 
hole Sgr A*. In this model, the variability of the line flux results 
from the propagation of an X-ray light front emitted by Sgr A* 
more than ~ 100 years ago. The discovery of an apparent super- 
luminal motion of the 6.4 keV line emission from the so-called 
"bridge" region provides strong support for this model (Ponti et 
al. 120101 ). The observed line flux variability with a timescale of 
a few years is hard to explain by a model of CR irradiation. 

In contrast to these results, the flux of the neutral or low- 
ionization Fe Ka line emitted from the Arches cluster vicinity 
does not show any significant variation over more than eight 
years of XMM-Newton repeated observations performed be- 
tween 2000 and 2009 (Sect. 4). Capelli et al. (1201 lbl l divided 



the zone of 6.4 keV line emission around the cluster into two 
subregions of about one parsec scale (labeled "N" and "S" by 
these authors) and found that both subregions emit the line at a 
constant flux. Other regions in the central molecular zone have 
been observed to emit a steady 6.4 keV line emission during 
about the same period, but they generally have larger spatial ex- 
tents (see, e.g., Ponti et al. [2010). Thus, the spatially averaged 
Fe Ka emission from Sgr B2 appears to have almost been con- 
stant for more than about seven years before fading away (Inui 
et al. 120091 Terrier et al. 120101 ). which is compatible with the 
light crossing time of the molecular cloud complex (see Odaka 
et al. 1201 11 1. But recent observations of Sgr B2 with Chandra 
suggest that the overall emission of the complex at 6.4 keV is in 
fact composed of small structures that have constantly changed 
shape over time (Terrier et al., in preparation). 

Together with the nondetection of time variability, the poor 
correlation of the spatial distribution of the 6.4 keV line emission 
with that of the molecular gas also argues against any origin of 
the Fe line in the Arches cluster region related to Sgr A* (Wang 
et al. 120061 . Lang et al. ( 120011 |2002 ) studied the position of 
the molecular clouds in the vicinity of the cluster by combining 
CS(2-1) observations with H92a recombination line data. The 
latter were used to trace the Arched filaments H n regions, which 
are thought to be located at edges of molecular clouds photoion- 
ized by the adjacent star cluster. Lang et al. (2001 , 2002) show 
that the molecular material in this region has a finger-like distri- 
bution and that the cluster is located in the midst of the so-called 
"-30 km s _1 cloud" complex, which extends over a region of 
~20 pc diameter (see also Serabyn & Guesten fl987l l. 

Figure Q3] compares the distribution of the 6.4 keV line 
emission around the star cluster with a high-resolution image 
in the hydrogen Paschen-a (Pa) line recently obtained with 
the Hubble Space Telescope/NICMOS instrument (Wang et al. 
120101 Dong et al. 1201 11 1. The Pa line emission is a sensitive 
tracer of massive stars - the Arches cluster is clearly visible in 
this figure at Galactic coordinates ((, b) ~ (0.122°, 0.018°) - and 
of warm interstellar gas photoionized by radiation from these 
stars. The main diffuse Pa-emitting features in Fig. [15] are the 
three easternmost Arched filaments: El (at 0.13° i £ i 0.15° 
and b ~ 0.025°), E2 (at b > 0.04°), and GO.10+0.02 (running 
from ((,b) ~ (0.09°, 0.006°) to (0.1°,0.025°); see also Lang 
et al. I20021 ). The 6.4 keV line emission is not well correlated 
with the Arched filaments. The most prominent structure at 6.4 
keV is concentrated in a region of only a few pc 2 surround- 
ing the star cluster, much smaller than the spatial extent of the 
-30 km s~' cloud complex. The origin of the faint Fe K emis- 
sion at ({, b) ~ (0.10°, 0.02°) is discusssed in the next section. 
This strongly suggests that the origin of the bright nonthermal 
X-ray radiation is related to the cluster itself and not to a distant 
source such as Sgr A*. 

Assuming that the nonthermal emission from the cloud re- 
gion is produced by a hard X-ray photoionization source located 
in the Arches cluster, the 4-12 keV source luminosity required 
to produce the observed 6.4 keV line flux can be estimated from 
Sunyaev & Churazov ([19981 : 
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where the distance to the GC is again assumed to be 8 kpc. 
Here, F 6.4 keV is the measured 6.4 keV line flux (Table [3), 
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is the column density of the line-emitting cloud, and Q the frac- 
tional solid angle that the cloud subtends at the X-ray source. 
This quantity is called the covering factor in, e.g., Yaqoob et al. 
(2010). In comparison, the unabsorbed luminosity of the cluster 
that we measured from the time-averaged XMM-Newton spec- 
tra is L x (4 - 12 keV) « 5 x 10 33 erg s _1 . Capelli et al. (12011 al l 
recently detected a 70% increase in the X-ray emission of the 
Arches cluster in March/ April 2007. However, the observed X- 
ray luminosity of the cluster is about two orders of magnitude 
short of what is required for the fluorescence interpretation. 

An alternative hypothesis is that the 6.4 keV line is produced 
by a transient photoionization source that was in a long-lasting 
(> 8.5 years) bright state at Lx ~ 10 36 erg s" 1 before a space tele- 
scope was able to detect it. No such source was detected with the 
Einstein observatory in 1979 (Watson et al. 1981) and with sub- 
sequent X-ray observatories as well, which imposes a minimum 
distance of ~ 4.6 pc between the cloud emitting at 6.4 keV and 
the putative transient X-ray source. This distance is increased to 
~ 9.2 pc if the cloud and the source are assumed to be at the 
same line-of-sight distance from the Earth. Furthermore, except 
for the extraordinarily long outburst of GRS 1915+105, which 
is predicted to last at least ~ 20 + 5 yr (Deegan et al. 2009|, the 
outburst duration of transient X-ray sources is generally much 
shorter than 8.5 years (see Degenaar et al. 2012). We also note 
that the Arches cluster is probably too young (f ~ 2.5 yr; Najarro 
et al. 120041) for an X-ray binary system to have formed within it. 

Thus, the 6.4 keV line emission arising from the vicinity 
of the Arches cluster is unlikely to result from photoionization 
and is most probably produced by CR impact. We have shown 
that the measured slope of the nonthermal power-law continuum 
(F = 1.6*?' 3 ,) and the EW of the 6.4 keV line from this region 
(EW6.4keV = 1.2+0.2 keV) are consistent with the predictions of 
the LECR ion model. On the other hand, LECR electrons cannot 
satisfactorily account for this emission, because it would require 
too high metallicity of the ambient gas (Z > 3.1 Z ) and too 
low minimum energy E m i n < 41 keV (Table @}. It is indeed un- 
likely that quasi-thermal electrons of such low energies can es- 
cape their acceleration region and penetrate a neutral or weakly 
ionized medium to produce the 6.4 keV line. We thus conclude 
that the 6.4 keV line emission from the cloud region is most 
likely produced by LECR ions. 

6.2.2. Are there other processes of production of the 6.4 keV 
line at work in the Arches cluster region? 

A relatively weak line at 6.4 keV is also detected in the spectrum 
of the X-ray emission from the star cluster. The low EW of this 
line (EW6.4 keV = 0.4+0. 1 keV) may suggest that this radiation is 
produced by LECR electrons accelerated within the cluster. The 
existence of a fast electron population there is supported by the 
detection with the VLA of diffuse nonthermal radio continuum 
emission (Yusef-Zadeh et al. 2003). The nonthermal electrons 
are thought to be produced by diffuse shock acceleration in col- 
liding wind shocks of the cluster flow. 

It is, however, more likely that the 6.4 keV line detected 
from this region is produced in molecular gas along the line of 
sight outside the star cluster. In the Chandra/ AOS 6.4 keV line 
image of the Arches cluster region, the bow shock-like struc- 
ture observed in the neutral or low-ionization Fe Ka line cov- 
ers the position of the star cluster (Wang et al. 20061. The 6.4 
keV line emission from the region "Cluster" is not observed in 
the present map (Figs. [10] and \15[ , because it has been artifi- 
cially removed in the process of subtraction of the continuum 
under the line (see Sect. 5.2). Lang et al. ( 120021 ) find evidence 




Fig. 15. XMM-Newton/EPIC continuum-subtracted 6.4-keV line 
intensity contours (linearly spaced between 3 x 10~ 8 and 
1.8 x 10~ 7 photons crrT 2 s" 1 arcmin" 2 ) overlaid with an 
aSr/NIC MOS m ap in the H Paschen-a line (Wang et al. l20T0l 
Dong et al. 1201 II) . The axes of the map indicate Galactic coordi- 
nates in degrees. The black ellipse and the magenta circle show 
the two regions used for spectral extraction (see Fig. [Tol l. The 
red arrow illustrates the observed proper motion of the Arches 
cluster, which is almost parallel to the Galactic plane (Stolte et 
al. 120081 Clarkson et al. 120121) . North is up and east to the left. 



of molecular gas lying just in front of the ionized gas associated 
with the most eastern Arched filament (El) close to the clus- 
ter sight line. According to the geometric arrangement of the 
-30 km s" 1 clouds proposed by these authors, it is likely that 
the cluster is presently interacting with this foreground molecu- 
lar gas. From the gradient of visual extinction detected by Stolte 
et al. (120021) over a field of 40" X 40" around the star cluster, 
9 < AAy < 15 mag, the H column density of this cloud along 
the line of sight can be estimated as i 3 x 10 22 crrT 2 . The 
calculations of the present paper show that LECR ions can pro- 
duce a significant 6.4 keV Fe Ka line emission in such a cloud 
(Sect. 2), especially in the case of strong particle diffusion for 
which the CR path length A can be much higher than (see 
Appendix A). It is thus likely that the weak nonthermal X-ray 
emission detected in the cluster spectrum has the same physical 
origin as the nonthermal emission from the cloud region. 

Relatively faint, diffuse emission in the neutral or low- 
ionization Fe Ka line is also detected to the west of the 
Arches cluster, from an extended region centered at ((, b) ~ 
(0.1°, 0.016°) (see Figs. [B]). Capelli et al. d201 lbb found the 
light curve of the 6.4 keV line flux from this region (labelled 
"SN" by these authors) to be constant over the 8-year observa- 
tion time. With a measured proper motion of ~ 4.5 mas yr" 1 al- 
most parallel to the Galactic plane and towards increasing longi- 
tude (Stolte et al. 120081 Clarkson et al. 120 121) . the Arches cluster 
was located within this region of the sky ~2x 10 4 years ago. It 
is therefore conceivable that this emission is also due to LECR 
ions that were accelerated within or close to the cluster at that 
time. That the nonthermal X-ray emission is still visible today 
would then indicate that the fast ions have propagated since then 
in a medium of mean density «h ~ 10 3 cm" 3 . Indeed most of the 
6.4 ke V line emission from LECR ions is produced by protons of 
kinetic energies < 200 MeV (see Fig. |5j)) and the slowing-down 
time of 200-MeV protons in «h = 10 3 cm" 3 is ~ 2 x 10 4 years. 
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Capelli et al. (2011b]) also considered a large region of 
6.4 keV line emission located at (t,b) ~ (0.11°, 0.075°) (see 
Fig- Q23 left panel). They measured a fast variability of the neu- 
tral or low-ionization Fe Ka line from this region and suggest 
that it could result from the illumination of a molecular cloud by 
a nearby transient X-ray source. The X-ray emission from this 
region is not studied in the present paper. 



7. Origin of the LECR ion population 

Two sites of particle acceleration in the Arches cluster region 
have been proposed. As already mentioned, Yusef-Zadeh et al. 
(2003 ) report evidence of diffuse nonthermal radio synchrotron 
emission from the cluster and suggest that the emitting relativis- 
tic electrons are accelerated by diffuse shock acceleration in the 
colliding stellar winds of the cluster flow. Another scenario is 
proposed by Wang et al. (120061) . who suggest that the 6.4 keV 
line emission from this region comes from LECR electrons pro- 
duced in a bow shock resulting from an ongoing supersonic col- 
lision between the star cluster and an adjacent molecular cloud. 
Both processes could also produce LECR ions. 

Since the work of Wang et al. ([2006), the apparent proper 
motion of the Arches cluster in the plane of the sky has been 
observed with Keck laser-guide star adaptive optics (Stolte et 
al. 12008 ; Clarkson et al. I2012I ). The direction of motion of the 
cluster stars relative to the field population is represented by the 
arrow in Fig. Q3] The 6.4 keV line emission close to the cluster 
shows two bright knots connected by a faint bridge to the east of 
the cluster, i.e. ahead of the moving stars (Fig. [TBI . The overall 
structure indeed suggests a bow shock. However, the Fe K line 
intensity scales as the product of the density of cosmic rays and 
that of the ambient medium around the cluster, which is probably 
highly inhomogeneous. A clear bow shock shape is therefore not 
to be expected. In fact, the 6.4 keV map from this region may 
also be explained by LECR ions escaping from the cluster and 
interacting with adjacent molecular gas. Thus, the morphology 
of the bright structure at 6.4 keV does not allow us to favor one 
of the two proposed sites for the production of fast ions. But 
more information can be obtained by studying the CR power 
required to explain the X-ray emission (Sect. 7.1), as well as the 
accelerated particle composition (Sect. 7.2). 



7.1. CR spectrum and energetics 

Whether the main source of LECR ions in the Arches cluster re- 
gion is the cluster bow shock or colliding stellar winds within the 
cluster flow, the nonthermal particles are likely to be produced 
by the diffusive shock acceleration (DSA) process. The nonther- 
mal particle energy distribution resulting from this process can 
be written for linear acceleration as (e.g. Jones & Ellison [T991l ) 



VJcs is the upstream sonic Mach number of the shock, whose 
velocity is V s . The sound velocity in the upstream gas is 
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Here, y g is the adiabatic index of the thermal gas upstream the 
shock front (y g = 5/3 for an ideal nonrelativistic gas) and Ms = 



where k is the Boltzmann constant, T the gas temperature and 
pmu the mean particle mass. In an interstellar molecular gas of 
temperature T = 100 K, c$ ~ 0.8 km s _1 . For a strong shock 
verifying V, » cs, we find from Eq. I© s^sa — 2, such that 
the particle spectrum in the nonrelativistic domain is a power 
law in kinetic energy of index s = 1.5 (see Eq. (0). Nonlinear 
effects due to the modification of the shock structure induced 
by the back-reaction of accelerated ions can slightly steepen the 
LECR spectrum, such that typically 1 .5 < s < 2 (Berezhko & 
Ellison IT999I ). The slope of the CR source spectrum that we de- 
rived from the X-ray spectral analysis, s = (see Tabled, 
is consistent with this theory. 

The total power acquired by LECR ions in the cloud region 
can be estimated from the best-fit normalization of the non- 
thermal X-ray component (Nlecr, see Table |4}. We find that 
the power injected by fast primary protons of energies between 
£Wn = 10 MeV and E max = 1 GeV in the X-ray emitting region 
is (4.3^2) x 10 38 erg s _1 (still assuming a distance to the GC 
of 8 kpc). Taking the uncertainty in E m [ n into account changes 
the proton power to (0.2-l)xl0 39 erg s _1 (see Sect. 5. 5). By in- 
tegration of the CR source spectrum, we find that about 30-60% 
more power is contained in suprathermal protons with E < E m i n , 
which, by assumption, do not penetrate dense regions of nonther- 
mal X-ray production. Considering the accelerated a-particles 
with C a /Cp = 0.1 adds another factor of 40%. The required to- 
tal CR power finally amounts to (0.5-1.8)xl0 39 erg s _1 . 



7.1 .1 . Mechanical power available from massive star winds 

The total mechanical power contained in the fast winds from 
massive stars of the cluster can be estimated from near in- 
frared and radio data. Using such observations, Rockefeller et 
al. ([2005) modeled the diffuse thermal X-ray emission from the 
cluster with 42 stellar wind sources with mass-loss rates in the 
range (0.3-17) X 10~ 5 M Q yr _1 and a terminal wind velocity of 
1000 km s _1 . The total mechanical power contained in these 42 
sources is 4 x 10 38 erg s -1 . Of course, only a fraction of this en- 
ergy reservoir can be converted to CR kinetic energy. We also 
note that LECR ions produced in the cluster are likely to diffuse 
away isotropically, such that those interacting with an adjacent 
molecular cloud emitting at 6.4 keV would probably represent a 
minority. Thus, the cluster wind is likely not powerful enough to 
explain the intensity of the nonthermal X-ray emission. 

7.1 .2. Mechanical power available from the Arches cluster 
proper motion 

The proper motion of the Arches cluster relative to the field star 
population has recently been measured to be 172 + 15 km s 
(Stolte et al. 2008 1 Clarkson et al. 2012 ). The cluster is also mov- 
ing away from the Sun, with a heliocentric line-of-sight veloc- 
ity of +95 + 8 km s (Figer et al. 120021 ). The resulting three- 
dimensional space velocity is V* « 196 km s _1 . To model the 
form of the bow shock resulting from this supersonic motion, 
we approximate the cluster as a point source object that loses 
mass at a rate Mw = 10~ 3 M Q yr~' through a wind of termi- 
nal velocity V w = 1000 km s _1 (see Rockefeller et al. [2005). 
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The shape of the bow shock is determined by the balance be- 
tween the ram pressure of the cluster wind and the ram pressure 
of the ongoing ISM gas. The pressure equilibrium is reached in 
the cluster direction of motion at the so-called standoff distance 
from the cluster (see, e.g., Wilkin 1996) 



Pbs - 



M W V W 



0.5 



= 2.4 pc , 



(6) 



where pic = 1 Am p n\c is the mass density of the local ISM. Here, 
we assume that since the birth of the cluster ~2.5 Myr ago (Figer 
et al. 120021 Najarro et al. 120041) . the bow shock has propagated 
most of the time in an intercloud medium of mean H density 
«ic ~ 10 cirT 3 (see Launhardt et al. l2002l for a description of the 
large-scale ISM in the GC region). 

The circular area of a bow shock projected on a plane per- 
pendicular to the direction of motion is At, s ~ !0nR^ s (see 
Wilkin 1 19961 1. Thus, the mechanical power processed by the 
cluster bow shock while propagating in the intercloud medium 
is Pic = 0.5picV 3 Ab s ~ 1.5 x 10 38 erg s _1 . In comparison, the 
steady state, mechanical power supplied by supernovae in the 
inner ~200 pc of the Galaxy is ~ 1.3 x 10 40 erg s _1 (Crocker et 
al. 1201 11 ). LECRs continuously accelerated out of the intercloud 
medium at the Arches cluster bow shock possibly contribute 
~1% of the steady-state CR power in the GC region (assuming 
the same acceleration efficiency as in supernova remnants). 

The initial total kinetic energy of the cluster motion is 
0.5M.V 2 ~ 1.9 x 10 52 erg, where M« ~ 5 x 10 4 M is the clus- 
ter initial total mass (Harfst et al. 2010). This energy would be 
dissipated in ~4 Myr according to our estimate of Pic- 

Most of the interstellar gas mass in the Galactic nuclear 
bulge is contained in dense molecular clouds with average H 
densities of «mc ~ 10 4 cirT 3 and a volume filling factor of 
a few percent (Launhardt et al. 2002). In the region where 
the Arches cluster is presently located, the volume filling fac- 
tor of dense molecular gas is even i 0.3 (Serabyn & Guesten 
I1987I ). Thus, the probability of a collision between the clus- 
ter bow shock and a molecular cloud is strong. The evidence 
that the cluster is presently interacting with a molecular cloud 
has already been discussed by Figer et al. (2002 ) and Wang et 
al. ((2006). This molecular cloud was identified as "Peak 2" in 
the CS map of Serabyn & Guesten (1987), who estimated its 
mass to be Mmc = (6 + 3) x 10 4 M and mean H density as 
«mc = (2 ± 1) x 10 4 cirT 3 . The corresponding diameter for a 
spherical cloud is d^c ~ 5.5 pc or 2.4' at a distance of 8 kpc, 
which is consistent with the apparent size of the cloud (Serabyn 
& Guesten [19"87l 

The total kinetic power processed in this collision is given 

by 



1 , 
3 mc = -Pmc(V. + Vmc) A c ■ 



(7) 



where pmc - l-4-m p riMC, Vmc * s tne velocity of the molecular 
cloud projected onto the direction of motion of the Arches clus- 
ter, and Ac the area of the contact surface between the "Peak 2" 
cloud and the bow shock. The latter quantity is not well known. 
We assume that it is equal to the area of the large region around 
the cluster emitting in the 6.4 keV line (i.e., the region labeled 
"Cloud" in Fig.[[0]and Tabled: A c = 7 pc 2 . 

The cloud-projected velocity Vmc obviously depends on the 
orbital path of the molecular cloud about the GC. By study- 
ing the velocity field of the molecular gas around the Arches 
cluster, Lang et al. ( 120011 12002) obtained constraints on the 
trajectory of the -30 km s _1 clouds. They find that the cloud 



complex likely resides on the far side of the GC, either on a 
X2 orbit (a noncircular orbit family set up in response to the 
Galaxy's stellar bar) or on a trajectory directed towards Sgr A* 
(the cloud complex would then be radially infalling into the su- 
permassive black hole) or perhaps on a trajectory midway be- 
tween the two situations. If the "Peak 2" cloud resides on an 
X2 orbit in the back of the Galaxy, the collision of this cloud 
with the Arches cluster is almost frontal (see Stolte et al. 2008 1 
Clarkson et al. 2012), and Vmc is close to the xj_ orbital speed 
v orb ~ 80 km s" 1 (see, e.g., Molinari et al. 1201 11 1. But if the 
cloud is radially infalling towards the supermassive black hole, 
given the radial velocity of the cloud (v, a d ~ -30 km s ) and the 
radial and transverse velocity components of the Arches clus- 
ter (v ra d =s +95 km s _1 and v tlans » +172 km s _1 directed to- 
wards positive longitude; see Clarkson et al. 2012 ), one finds that 
Vmc ~ 20 km s _1 . Thus, depending on the exact cloud trajectory 
Vmc ~ 50 ± 30 km s -1 , which, together with V, as 196 km s , 
gives P M c ~ 2.3 x 10 40 erg s _1 from Eq. (0. 

In comparison, the CR power needed to explain the X-ray 
observations is dW/dt = (0.5-1.8)xl0 39 erg s , such that the 
required particle acceleration efficiency in the bow shock sys- 
tem amounts to a few percent. This is a typical efficiency in 
the DSA theory and in the phenomenology of the acceleration 
of the Galactic CRs in supernova remnant shocks, as well (see, 
e.g., Tatischef fl2008l and references therein). However, a detailed 
study of the particle acceleration process at work in this peculiar 
shock system would go beyond the scope of this paper. 

The collision kinetic power estimated above is comparable 
to the steady-state, mechanical power due to supernovae in the 
inner Galaxy, ~ 1.3 X 10 40 erg S -1 (Crocker et al. I20TT1 ). But 
the typical time duration of a collision between the Arches clus- 
ter bow shock and a molecular cloud is expected to be only 
~ 3 x 10 4 yr, assuming that the size Jmc ~ 5.5 pc is typical of the 
highly-fragmented dense molecular gas of the GC region. The 
mechanical energy released in such a collision is then ~10 52 erg, 
i.e. comparable to the cluster initial total kinetic energy. Thus, 
the cluster bow shock very likely has collided no more than 
once with a molecular cloud since the cluster birth ~2.5 Myr 
ago. Such a collision can briefly release in the ISM a power in 
LECRs comparable to the steady-state CR power supplied by 
supernovae in the GC region (Crocker et al. 1201 11 1. 

The ISM volume swept up by the Arches cluster bow shock 
since the cluster birth can be estimated as 



V bs = A bs V«f* 



9 x 10 4 pc 3 



(8) 



where f„ =2.5 Myr is the estimated cluster age. This volume rep- 
resents less than 1 % of the total volume of the Galactic nuclear 
bulge, Vnb ~ 1.5 x 10 7 pc 3 . In comparison, the volume filling fac- 
tor of dense molecular cloud in the inner ~230 pc of the Galaxy 
is a few percent (Launhardt et al. 2002). It is thus likely that the 
star cluster did not experience any interaction with a molecular 
cloud before the one with the "Peak 2" cloud that is presently 
observed. 

Simulated orbits of the Arches cluster about the GC suggest 
that the cluster formed in the front of the Galaxy near an xi orbit 
(Stolte et al. |2008 ). That the Arches cluster presently interacts 
with a molecular cloud located behind the GC shows that the 
cluster's orbit is retrograde to the general motion of stars and 
gas clouds in the bar potential (see Fig. 8 of Stolte et al. 2008). 
According to the simulation of possible orbits, the cluster has 
performed about half a revolution around the GC since its for- 
mation, which may have brought it near the far side of the el- 
liptical ring of dense molecular clouds recently studied with the 
Herschel satellite (Molinari et al. 1201 11 1. In this environment, the 
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probability of a collision between the cluster and a molecular 
cloud has become strong. 

7.2. Accelerated ion composition 

Fast C and heavier ions can emit very broad X-ray lines result- 
ing from 2p to Is (Ka) and 3p to Is (K/3) in-flight transitions. 
The 2p and 3p orbital states can be populated either by electron 
capture from ambient atoms (i.e. charge exchange) or by excita- 
tion of Is electrons for fast ions having one or two electrons. To 
study the composition of the energetic ions accelerated near the 
Arches cluster, we developed new LECR ion models that include 
the line emission of fast C, N, O, Ne, Mg, Si, S, and Fe. We used 
the tables of K X-ray differential multiplicities, dMf k /dE, given 
in Tatischeff et al. ( 1998). This quantity is defined as the number 
of photons emitted in the Kk line by the projectile i as it slows 
down over the differential kinetic energy interval dE, owing to 
interactions with all the constituents of the ambient medium. In 
the adopted steady-state, slab interaction model, the X-ray line 
production rate is then simply given by 

dQf k f 00 dMf k f £ A< £ > do, 

-%-(Ex)= —±r(E x ,E)dE -^(E')dE'. (9) 

dt Jo dE J E dt 

The energy of the emitted X-rays depends on both the veloc- 
ity and spatial distributions of the fast ions through the usual 
Doppler formula. We assumed isotropic propagation of the 
LECRs in the interaction region, which leads to a maximum 
broadening of the lines. 

We note that the calculations of Tatischeff et al. ( 1998) were 
done for an ambient medium of solar composition. Nevertheless, 
their multiplicity results can be used in good approximation for 
a medium of order twice solar metallicity, since it was found 
that most of the line emission from the fast ions is produced by 
interactions with ambient H and He. 

We first developed an X-ray production model in which the 
abundances of the heavy ions are in solar proportions relative to 
each other, but can vary with respect to H and He, that is (see 
also Eq. SKJ\ ) 

^=fj£r) ■ (10) 

The parameter / met adds to the spectral index s and the metallic- 
ity of the ambient medium Z as a free parameter of the model. 
But here, we fixed E m m = 10 MeV and A = 5 x 10 24 cm~ 2 , as 
in Sect. 5.2. The resulting model was made readable in XSPEC 
and then used to fit the XMM-Newton spectra extracted from the 
cloud region. As before, we chose to fit to the data the global 
XSPEC model WABSx(APEC + LECR/). The best fit was ob- 
tained for / met = 3.9*i A The corresponding^ 2 and best-fit val- 
ues of the other parameters are nearly the same as in Table |U 
in particular with Z/Z Q = 1.7 ± 0.2 and s = 2.0+° *. The ob- 
tained limit at the 90% confidence level / met < 15.2 shows that 
the heavy ion abundance is not well constrained by this model. 

Figure [T6h shows calculated X-ray spectra for the solar sys- 
tem composition with / met = 1 and 10. We see that intense 
broad lines can be emitted below 1 keV from X-ray transi- 
tions in fast C, N, and O. But this line emission cannot be ob- 
served from sources in the GC region, because of the strong in- 
terstellar photoelectric absorption along the line of sight. With 
Nu ~ 10 23 H cm~ 2 (Table |U), the main constraint on the ac- 
celerated particle composition is provided by a very broad line 
feature between ~5.5 and 9 keV owing to de-excitations in fast 
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Fig. 16. Calculated X-ray emission produced by LECR ions with 
the spectral parameters s = 1.9, E m [ n = 10 MeV nucleon" 1 , 
and the escape path length A = 5 x 10 24 crrT 2 , interacting in a 
gas cloud of metallicity Z = 1.7Z (as obtained in Sect. 5). In 
panel a) the abundances of the accelerated heavy ions are in solar 
proportion, but with a metallicity enhancement factor relative to 
the solar system (SS) composition (see text) / met = 1 (dashed 
curve) and f met = 10 (solid and dotted curves). Panel b): same 
but for the CRS composition (see text). The dotted curves show 
the effect of photoelectric absorption on the X-ray spectra for a 
H column density of 10 23 cm 2 . 



Fe. This emission is produced by Fe ions of energies between 
~5 and 20 MeV nucleon 1 (Tatischeff et al. [19981 . 

In a second model, we assumed that the fast metals have the 
composition as the current epoch Galactic CRs at their sources. 
We obtained the CR source (CRS) composition by taking the 
heavy ion abundances relative to O given by Engelmann et 
al. (T990) and using the abundance ratios C a /Co = 19 and 
C p /C a = 15 recommended by Meyer et al. ( 1997). The resulting 
CRS composition is consistent with the recent theoretical works 
of Putze et al. d201 1) . The Fe abundance in the CRS composition 
is C-pe/Cp = 6.72 x 10~ 4 , which is 19 times higher than the one in 
the solar system composition, (Cp e /C P ) G = 3.45 x 10~ 5 (Lodders 
et al. 120031 . The best fit of this model to the X-ray spectra of the 
cloud region was obtained for / met = O.ll^ f? (i.e. / met < 0.71 at 
the 90% confidence level), consistent with the higher abundance 
of Fe in the CRS composition. X-ray spectra calculated for this 
composition are shown in Fig.[T6b. 

The main outputs of this analysis are the nondetection in the 
X-ray spectra of the cloud region of a significant excess emis- 
sion from fast Fe and the implication that Cp e /C p i 5 x 10~ 4 . 
This result by itself does not provide strong support for one 
or the other possible site of acceleration of the LECR ions in 
the Arches cluster region. Indeed, with the best-fit metallicity 
Z/Z = 1 .7 + 0.2, the Fe abundance in the local molecular cloud 
is aFe = Z/Z (flFe)o ~ 5.9 x 10~ 5 , which is well below the upper 
limit obtained above. The Fe abundance can be slightly higher 
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in the wind material expelled by the massive stars of the Arches 
cluster. According to Parizot et al. (1997), one expects in the 
average composition of the winds from OB associations in the 
inner Galaxy ape = 1.4 x 10~ 4 , which is still below the derived 
upper limit. 

However, the Galactic CRS composition is best described 
in terms of a general enhancement of the refractory elements 
such as Fe relative to the volatile ones (Meyer et al. 1997 ). This 
selection effect is most likely related to the acceleration pro- 
cess at work in supernova remnant shock waves. Given the limit 
/met < 0.71 that we obtained for the CRS composition, we con- 
clude that this effect is weaker in the shock system associated 
with the Arches cluster proper motion, which according to the 
energetics arguments discussed in Sect. 6.2.2, is the most likely 
site of acceleration of the X-ray emitting LECR ions. 

8. Cosmic-ray ionization rate 

We estimated in Sect. 7.1 that a kinetic power of (0.5- 
1.8)xl0 39 erg s _1 is currently delivered by the Arches cluster 
bow shock system to LECR ions of energies < 1 GeV nucleon -1 . 
The power continuously deposited into the adjacent molecular 
cloud is lower than that, because of (i) the nonpenetration of 
CRs with E < E m \ n into the interaction region and (ii) the escape 
from the cloud of the highest energy CRs. For A = 5 x 10 24 cm" 2 , 
protons of energies up to 180 MeV are stopped in the cloud, 
whereas those injected at higher energies do not virtually lose 
energy in this medium. Taking these effects into account, we find 
that the power deposited by LECRs into the cloud amounts to 
Wd ~ 4 x 10 38 erg s _1 . The initial kinetic energy of the fast ions 
is essentially lost through ionization of the ambient gas and the 
corresponding ionization rate can be estimated to be 

l.4m p Wd , , 

<Th = — ~ — ~ 10 H s , (11) 

eiM M c 

where e, » 40 eV is the mean energy required for a fast ion 
to produce a free electron in a neutral gas mixture of H2 and 
He in cosmic proportion (Dalgarno et al. 1999) and Mmc = 
(6 + 3) x 10 4 M s is the cloud mass (Serabyn & Giiesten 1987 ). 
The mean ionization rate induced by LECRs in this molecular 
cloud is significantly higher than the mean ionization rate in the 
GC region, £ H ~ 10~ 15 H 1 s" 1 (see Crocker et al. l20TT1 and ref- 
erences therein). 

By integrating the differential equilibrium number of LECRs 
in the X-ray production region (see Eq. ( IA.U ). we find that the 
total kinetic energy contained in fast ions diffusing in the cloud 
is E tot ~ 4 x 10 48 erg. The corresponding mean energy density 
is E tQ[ /V M c ~ 1000 eV cirT 3 (here V MC = M MC /(lAn MC m p )), 
which is about one thousand times higher than the Galactic CR 
energy density in the solar neighborhood. Thus, the molecular 
cloud irradiated by fast particles accelerated near the Arches 
cluster bow shock shows some similarities with the "extreme 
CR dominated regions" recently studied by Papadopoulos et al. 
( 1201 11 1 in the context of starbursts. Following the works of these 
authors, LECRs could explain the high temperature of the "Peak 
2" cloud measured by Serabyn & Giiesten JT987T ): T > 100 K. 

9. Gamma-ray counterparts 

Collisions of LECR ions with molecular cloud matter can lead to 
nuclear excitations of both ambient and accelerated heavy ions 
followed by emission of de-excitation y-ray lines (Ramaty et al. 



[19791 Benhabiles-Mezhoud et al. 1201 11 1. We calculated the y- 
ray line flux expected from the Arches cluster region using the 
same CR interaction model as before (see Appendix A), with 
s = 1.9, E min = 10 MeV, JVlecr = 5.6 x 10~ 8 erg cirT 2 s" 1 , 
A = 5 x 10 24 cm" 2 , and Z/Z = 1.7 (Table©. The predicted flux 
is well below the sensitivity limit of the INTEGRAL observatory. 
For example, we obtain a flux of 2.3 x 10~ 8 ph cm -2 s _1 in the 
4.44 MeV line from de-excitations of ambient 12 C, whereas the 
sensitivity of the INTEGRAL spectrometer SPI for detection of a 
narrow line at this energy is > 10~ 5 ph cirT 2 s _1 . 

Nuclear interactions of CR ions with ambient matter can also 
lead to high-energy y-ray emission via the production and subse- 
quent decay of jfi mesons. We used the model of Dermer ( 1986 ) 
for this calculation, but multiplied the 7r° emissivity given by this 
author by a factor of 1 .27 to be consistent with the local emissiv- 
ity measured with the Fermi Gamma-ray Space Telescope (Abdo 
et al. I2Q091 ). The high-energy y-ray flux strongly depends on the 
shape of the CR energy distribution, because the neutral pions 
are produced at significantly higher energies than the nonther- 
mal X-rays. We first assumed a CR source spectrum of the form 
given by Eq. (0 for this calculation (i.e. resulting from the DSA 
process) with sdsa = 2s — 1 = 2.8 and no high-energy cutoff. 
We then found that the Arches cluster region would emit a flux 
of 5.7 x 10~ 7 ph cirT 2 s in y-rays of energies > 300 MeV. 
Such high-energy emission would have probably been already 
detected by Fermi, since the predicted flux is ~ 1 .75 times higher 
than the flux of the Galact ic central source 1FGL J1745. 6-2900 
(Chernyakova et al. 1201 11 1. But an exponential cutoff in the CR 
distribution can be expected either because of the finite size 
of the particle acceleration region near the cluster bow shock 
or the finite time available for particle acceleration. For ex- 
ample, with an exponential cutoff at 0.5 GeV nucleon 1 (resp. 
1 GeV nucleon 1 ), the flux of y-rays > 300 MeV would be re- 
duced to 1.4 x 10~ 8 ph cirT 2 s _1 (resp. 5.1 x 10~ 8 ph cirT 2 s _1 ) 
without significantly changing the nonthermal X-ray production. 
The high-energy y-ray emission from the Arches cluster region 
would then be undetectable with Fermi. 

10. Summary 

We have studied the production of nonthermal line and contin- 
uum X-rays by interaction of LECR electrons and ions with a 
neutral ambient medium in detail. We developed a steady-state, 
slab model in which accelerated particles penetrate at a constant 
rate a cloud of neutral gas, where they produce nonthermal X- 
rays by atomic collisions until they either stop or escape from 
the cloud. We examined the properties of the neutral Fe Ka line 
excited by impacts of LECR electrons and ions. The predicted 
line EW and luminosity, as well as the slope of the underlying 
bremsstrahlung continuum, were presented as functions of the 
free parameters of the model. These results are intended to help 
observers study the potential role of LECRs for any 6.4 keV line 
emission and possibly decipher the nature of the nonthermal par- 
ticles responsible for the line emission. In addition, we generated 
LECR electron and ion models that can be used in the XSPEC 
software for more quantitative comparison with data. 

We showed, in particular, that the EW of the neutral Fe 
Ka line excited by LECR electrons is generally expected to 
be lower than 1 keV, except if the metallicity of the ambient 
medium exceeds *2Z . But LECR ions with a relatively soft 
source spectrum can lead to a much larger EW. However, the 
production of 6.4 keV line photons by both LECR electrons 
and ions is relatively inefficient: the radiation yield R(,a = 
L X (6A keV)/(dW/dt) is typically on the order of 10~ 6 , mean- 
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ing that a high power in LECRs should generally be needed to 
produce an observable neutral Fe Ka line. 

We then employed the newly developed models to study the 
X-ray emission emanating from the Arches cluster region. We 
used all public XMM-Newton EPIC observations encompassing 
the studied region for our analysis. The main results of this anal- 
ysis can be summarized as follows. 

- The X-ray flux detected from the Arches cluster is domi- 
nated by the emission of an optically thin thermal plasma 
with a temperature kT ~ 1.7 keV. This component most 
likely arises from the thermalization of massive star winds 
that merge and expand together, plus the contribution of sev- 
eral colliding stellar wind binaries within the cluster. 

- A second thermal plasma of lower temperature is required to 
explain the presence of He-like Si and S Ka emission lines 
in the X-ray spectrum of the cluster. This component, which 
was not detected in previous X-ray observations, may be pro- 
duced by a collection of individual massive stars. 

- Bright 6.4 keV Fe Ka line structures are observed around the 
Arches cluster. We found that the line flux from this region is 
consistent with its being constant over more than eight years 
of XMM-Newton repeated observations, in agreement with 
the recent works of Capelli et al. (2011b ). This radiation is 
unlikely to result from the photoionization of a molecular 
cloud by a hard X-ray source. It is also probably not pro- 
duced by LECR electrons, because it would require a metal- 
licity of the ambient gas (Z > 3.1 Z ) that is too high. On the 
other hand, the X-ray emission observed around the cluster 
can be well-fitted with a model composed of an optically thin 
thermal plasma and a nonthermal component produced by 
LECR ions. The best-fit metallicity of the ambient medium 
found with this model is Z/Z Q = 1 .7 ± 0.2, and the best-fit 
CR source spectral index is s = 1.9*2 5. 

- The required flux of LECR ions is likely to be produced by 
the diffusive shock acceleration process in the region of in- 
teraction of the Arches cluster and the adjacent molecular 
cloud identified as "Peak 2" in the CS map of Serabyn & 
Giiesten (1987). We estimated that a total kinetic power of 
~ 2.3 x 10 40 erg s -1 is currently processed in the ongoing 
supersonic collision between the star cluster and the molec- 
ular cloud emitting the 6.4 keV line. A particle acceleration 
efficiency of a few percent in the resulting bow shock system 
would produce enough CR power to explain the luminosity 
of the nonthermal X-ray emission. 

- We developed LECR ion models that include the production 
of broad X-ray lines from fast C and heavier ions follow- 
ing electron captures from ambient atoms (i.e. charge ex- 
changes) and atomic excitations. It allowed us to constrain 
the abundance of fast Fe ions relative to protons in the LECR 
ion population: Cp e /C p i5 x 10~ 4 . This limit is ~15 times 
higher than the Fe abundance in the solar system composi- 
tion. 

- The mean ionization rate induced by LECRs in the molecular 
cloud that is thought to presently interact with the Arches 
cluster is £h ~ 10~ 13 H 1 s -1 . The CR energy density in 
the interaction region is estimated to be ~ 1000 eV cm -3 , 
which is about one thousand times higher than the Galactic 
CR energy density in the solar neighborhood. 

- The high-energy y-ray emission produced by hadronic col- 
lisions between CRs accelerated in the Arches cluster bow 
shock system and ambient material might be detected with 
the Fermi Gamma-ray Space Telescope. It crucially depends, 



however, on the unknown shape of the CR energy distribu- 
tion above ~1 GeV nucleoli -1 . 

The nonthermal X-ray emission emanating from the Arches 
cluster region probably offers the best available signature cur- 
rently for a source of low-energy hadronic cosmic rays in the 
Galaxy. Deeper observations of this region with X-ray telescopes 
would allow better characterization of the acceleration process 
and the effects of LECRs on the interstellar medium. The theory 
presented in this paper could also be useful for identifying new 
sources of LECRs in the Galaxy. 
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Appendix A: Cosmic-ray interaction model 

We consider a model in which low-energy cosmic rays (LECRs) 
are produced in an unspecified acceleration region and pen- 
etrate a nearby cloud of neutral gas at a constant rate (see 
Fig. IA. lb . The energetic particles can produce nonthermal X- 
rays by atomic collisions while they slow down by ionization 
and radiative energy losses in the dense cloud. We further as- 
sume that the LECRs that penetrate the cloud can escape from it 
after an energy-independent path length A, which is a free pa- 
rameter of the model. The differential equilibrium number of 
primary CRs of type i (electrons, protons, or a particles) in the 
cloud is then given by 



(A.l) 



i r^> dQi , , 

ME) = —(E')dE' , 

[dE/dt(E)] t J E dt 

where (dQijdt) is the differential rate of LECRs injected in the 
cloud, [dE /dt(E)]i is the CR energy loss rate, and the maximum 
energy E' A (E) is related to the escape path length A (expressed 
in units of H atoms crrT 2 ) by 



dE' 



[dE/dN K (E')]i ' 



where 
dE 



dN H j i 



1 (dE 
— I — m„ 

V;H H \ dt ) i 



dE\ IdE 
dx ) i H \ dx jj 



(A.2) 



(A.3) 



Here, v, is the particle velocity, «h the mean number density 
of H atoms in the cloud, m p the proton mass, one = 0.0964 
the cosmic abundance of He relative to H (Lodders 2003), and 
(dE/dx)i,H and (dE/dx)i t ue the CR stopping powers (in units of 
MeV g -1 cm 2 ) in ambient H and He, respectively. We used for 
electrons the stopping-power tables of Berger & Seltzer (1982) 
below 1 GeV and the relativistic formulae given by Schlickeiser 
( 2002 ) above this energy. The stopping powers for protons and 
ff-particles were extracted from the online databases PSTAR and 
ASTAR, respectively (Berger et al. [2005l . 

The process of CR transport in the cloud, which does not 
need to be specified in the above formalism, is nevertheless rel- 
evant to estimate the escape path length A from the cloud size. 
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It is clear that if the cloud medium is not diffusive, because of, 
e.g., efficient ion-neutral damping of MHD waves, A ~ n^Lc, 
where Lq is the characteristic size of the cloud. But otherwise, 
the escape path length, which can then be estimated as 



A 



L c Vin H 
6D 



(A.4) 



can be much greater than the characteristic column density iV£ = 
rinLc, depending on the diffusion coefficient D. For example, 
with the typical mean diffusion coefficient for the propagation 
of Galactic CR nuclei in the local interstellar magnetic field B 
(Berezinsky et al l 19901 ), 



D * 10 28 jS( 



B 



1GV/ \3yuG 



-0.5 



2 -1 

cm s , 



(A.5) 



where B = v,/c and R, is the particle rigidity, one gets from 
Eq. ( 1A.4) for non-relativistic protons: 



A ~ 5 x 10 24 



-0.25 



AfC 



50 MeV / 
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100 nG 
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(A.6) 



where A^, «h and B are scaled to typical values for massive 
molecular clouds in the GC region. 

For nonrelativistic particles diffusing in the cloud with a dif- 
fusion coefficient D oc BR*" typically with 0.3 < sq < 0.5, the 
escape path length estimated from Eq. ( 1A.4) depends only mildly 
on energy as A oc E^ SD ^ 2 . However, we have adopted here a sim- 
ple slab model with an energy-independent escape path length in 
order to limit the number of free parameters as much as possible. 

The process of CR penetration into molecular clouds is not 
well known (see, e.g., Gabici et al. 2007 and references therein). 
The theoretical predictions range from almost-free penetration 
(e.g. Cesarsky & Volk 1978 ) to exclusion of CRs of kinetic en- 
ergies up to tens of GeV (e.g. Skilling & Strong 1976). Here, for 
simplicity, we assume that CRs can freely penetrate the clouds if 
their kinetic energy is higher than a threshold energy E m i n , which 
is another free parameter of the model. We further consider the 
differential rate of primary CRs that penetrate the nonthermal X- 
ray production region to be a power law in kinetic energy above 



dQi_ 
dt 



(E) = C t E~ s for E > E n 



(A.7) 



The model finally has four free parameters that can be stud- 
ied from spectral fitting of X-ray data (see Sect. 5): A, E m [ n , 
the power-law spectral index s, and the metallicity of the X-ray 
emitting cloud, Z. The X-ray spectral analysis also provides the 
CR spectrum normalization C,-, which allows one to estimate the 
power injected by the primary LECRs into the nonthermal X-ray 
production region: 



dWj 
dt 



r 



E'—(E')dE' . 
dt 



(A.8) 



In the following, the integration in the above equation is limited 
to E m . dx = 1 GeV. Due to CR escape, the power continuously 
deposited by the fast particles inside the cloud should generally 
be lower than dWi/dt. 



Appendix B: X-rays from accelerated electron 
interactions 

In the framework of the adopted steady-state, slab interaction 
model, the differential X-ray production rate from collisions of 
accelerated electrons with the cloud constituents can be written 

as 



dt 



(E x ) = « H 







dcr. 



<-'.) 



dEy 



(E x ,E)v e (E) 



x [N e JE) + N e ,(E)]dE. 



(B.l) 



where aj is the abundance of element j relative to H in the X-ray 
emitting cloud, (do- e j/dEx) is the differential X-ray production 
cross section for electron interaction with atoms j, and N e4 , and 
N etS are the differential equilibrium numbers of primary and sec- 
ondary LECR electrons in the ambient medium, respectively. 

B. 1 . Secondary electron production 

Primary LECR electrons injected into an interstellar molecular 
cloud produce secondary electrons mainly from ionization of 
ambient H2 molecules and He atoms. The corresponding differ- 
ential production rate of knock-on electrons is given by 



dQ e , 
dt 



(E s ) 



«H 



r*O0 

JlE, 



0.5^^(E S ,E P ) +a He ——(E s ,E p ) 
dE, 



dE s 



x v e (E p )N e , p (E p )dE p 



(B.2) 



where (d<Tn 2 /dE s ) and (do-n e /dE s ) are the H2 and He differen- 
tial ionization cross sections for the production of a secondary 
electron of energy E s by impact of a primary electron of energy 
E p . The lower limit of the integral is 2E S , because the primary 
electron is by convention the faster of the two electrons emerg- 
ing from the collision. The maximum possible energy transfer is 
therefore E s = 0.5(£,, - Bj) =* 0.5E P , where B H2 = 15.43 eV and 
Z?He = 24.59 eV are the electron binding energies of H2 and He, 
respectively. This convention is consistent with the definition of 
the stopping powers used throughout this paper (see Eq. ( 1A.31 )). 
which also pertain to the outgoing electron of higher energy. 

The differential ionization cross sections are calculated from 
the relativistic binary encounter dipole (RBED) theory (Kim et 
al. 119941 12000bl ). This successful model combines the binary- 
encounter theory for hard collisions with the dipole interaction 
of the Bethe theory for fast incident electrons. For the differ- 
ential oscillator strengths, we use the analytic fits provided by 
Kim et al. (fT994l ) for H 2 and Kim et al. (12000a) for He. For the 
average orbital kinetic energy of the target electrons, we take 
Un 2 = 15.98 eV and U He = 39.51 eV. 

By inserting Eq. (IA.1) into Eq. ( IB. 2) and using for the elec- 
tron energy loss rate the expression given in Eq. (IA.3) , we see 
that the secondary electron production rate does not depend on 
the absolute density of H atoms in the ambient medium («h). 
This comment also applies to the X-ray production rate, which 
only depends on the relative abundances aj (see Eq. (IB.D ). 
This important property of the adopted steady-state, slab model 
will allow us to estimate unambiguously the cosmic -ray power 
dWi/dt (Eq. lA.8) from the measured X-ray flux. 

Calculated differential production rates of primary and 
knock-on electrons are shown in Fig. IB. II Also shown is the cor- 
responding steady-state differential number of secondary elec- 
trons in the ambient medium, N ClS . We calculated the latter from 
Eqs. (IA.1) and (IA.2) . assuming the characteristic escape path 
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Fig. B.l. Calculated differential equilibrium electron numbers 
(N e ; solid lines) for two differential injection rates of primary 
electrons (dQ e . p /dt; dotted lines): a) s-3, Zimi^lO keV; b) s=2, 
£min=l MeV. Also shown are the differential production rates 
of secondary, knock-on electrons (dQ^Jdt; dashed lines). The 
H density in the nonthermal X-ray production region, which in- 
tervenes in the calculation of N e , p and N eyS , is «h = 10 4 cm -3 
and the path length of the primary electrons in this region is 
A = 10 24 cm~ 2 . The calculations are normalized to a total power 
of 1 erg s injected by the primary LECR electrons in the X-ray 
production region. 



length of the secondary particles to be A/2. Although this as- 
sumption is uncertain, it has no significant effect on the total 
X-ray production. 

We see in Fig. IB. II that the effect of H2 and He ionization 
on the electron energy distribution is to redistribute the total ki- 
netic energy of the injected particles to a larger number of lower- 
energy electrons. Thus, for hard enough primary electron spec- 
trum (i.e. low s and high E m i n , see Fig. IB. lb ), secondary elec- 
trons of energies E s i 10 keV could potentially make a signifi- 
cant contribution to the total nonthermal X-ray emission. On the 
other hand, one can easily check that the successive production 
of knock-on electrons by the secondary electrons themselves can 
be safely neglected for the X-ray emission. 



B.2. X-ray continuum emission 

The X-ray continuum emission is due to the bremsstrahlung 
of both primary and secondary electrons. We take electron 
bremsstrahlung into account only in ambient H and He and cal- 
culate the differential cross sections from the work of Strong 
et al. (120001 appendix A), which is largely based on Koch & 
Motz (I19591 I. We use the scattering functions from Blumenthal 
& Gould (119701 1 to take into account the arbitrary screening of 
the H and He nuclei by the bound electrons. 
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Fig. B.2. Calculated X-ray emission produced by LECR elec- 
trons with the source spectra shown in Fig. IB. II interact- 
ing in a gas cloud of solar metallicity. PEB: primary elec- 
tron bremsstrahlung; SEB: secondary electron bremsstrahlung. 
Photoelectric absorption is taken into account with a H column 
density of 10 2 



-|22 cm 2 



B.3. X-ray line emission 

The X-ray line emission results from the filling of inner-shell 
vacancies produced by fast electrons in ambient atoms. We con- 
sider the Ka and K/3 lines (2p— >ls and 3/?— >ls transitions in the 
Siegbahn notation) from ambient C, N, O, Ne, Mg, Si, S, Ar, Ca, 
Fe, and Ni. The corresponding cross sections can be written as 



dcr Kk 

(T K HE) = -^(E X ,E) = 6(E X - E Kk )cr' ej (E)u>f k , 



dE 



(B.3) 



x 



where E X k is the energy of line K£ (Ka or K/3), 5(Ex - E^) is 
Dirac's delta function, cr^JiE) the total cross section for the K- 

shell ionization of atom j by an electron of energy E, a nd a>^ k 
the KA: fluorescence yield for atom j (Kaastra & Me we [19931 . 
Note that d K ^-Q for element /' with atomic number < 12 (i.e. 
Mg), since these atoms do not have 3p electrons in their ground 
level. 

For the K-shell ionization cross sections, we adopted the 
semi-empirical formula of Quarles (1976), which agrees well 
with the RBED cross sections for Ni and lighter elements (see 
Kim et al. I2000bb and is simpler to use. We checked that the 
Quarles's formula correctly reproduces the data compiled in 
Long et al. (1990), in particular at relativistic energies. 

The width of the X-ray lines produced by electron impact 
can be estimated from the sum of the natural widths of the 
atomic levels involved in the transition. Indeed, broadening ef- 
fects caused by multiple simultaneous ionizations can be safely 
neglected for LECR electrons. Thus, the Kai and Kct2 compo- 
nents of the Fe Ka line have experimental full widths at half- 
maximum (FWHM) of only 2.5 and 3.2 eV, respectively (Salem 
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& Lee 119761) . However, the energy separation of the two fine- 
structure components is 13 eV, which is much less than the en- 
ergy resolution at 6.4 keV of the X-ray cameras aboard XMM- 
Newton and Chandra, but larger than the expected resolution of 
the ASTRO-H X-ray Calorimeter Spectrometer (7 eV FWHM; 
Takahashi et al. 12010) . Here, we neglect the fine-structure split- 
ting of the K lines and for simplicity adopt the same width for 
all the lines: AE X = 10 eV. 

Figure lBT2l shows calculated nonthermal X-ray spectra (Lx = 
ExXdQx/dt) produced by LECR electrons injected with the dif- 
ferential rates shown in Fig. lB.ll into a cloud of solar metallicity. 
We took the photoelectric absorption of X-rays into account us- 
ing a H column density Nu = 10 22 cm" 2 and the cross sections 
of Morrison & McCammon ( I1983I ). We see in Fig. IB. 21 that the 
most prominent line is that of Fe at 6.40 keV. This is because 
this element has the highest product of Ka fluorescence yield 
(Wp" =0.3039, Kaastra & Mewe 1993 ) and cosmic abundance. 
The EW of the Fe Ka line is equal to 293 and 394 eV in the spec- 
tra shown in panels a and b, respectively. The second strongest 
line in these spectra is the Si Ka line at 1 .74 keV; its EW is equal 
to 80 and 90 eV in panels a and b, respectively. We also see in 
this figure that (i) the shape of the continuum emission reflects 
the hardness of the primary electron injection spectrum, and (ii) 
the total X-ray emission is dominated by the contribution of the 
primary electrons. The emission from the secondary electrons 
is negligible in panel a and accounts for 10-20 % of the total 
emission below 10 keV in panel b. 

Appendix C: X-rays from accelerated ion 
interactions 

The differential X-ray production rate from accelerated ion in- 
teractions can be written with a slight modification of Eq. ( IB. I) , 
as follows: 



10 



dQx , 
dt 



(E x ) = « H J] aj j J] 



+ ^r(Ex,E)v e (E)N e JE) 
dE x 



dE 



(CI) 



where the index i runs over the constituents of the nonthermal 
ion population. The first term in the integral represents the X- 
ray production by the primary LECR ions and the second term 
the contribution of the secondary electrons. As a starting point, 
we assume in the present work that the LECR ion population is 
mainly composed of protons and a particles and that the con- 
tributions of accelerated metals to the total X-ray emission can 
be neglected. We therefore do not consider the broad X-ray 
line emission that can arise from atomic transitions in fast C 
and heavier species following electron captures and excitations 
(Tatischeff et al. 119981 ), except in Sect. 7.2. However, for typical 
compositions of accelerated cosmic particles, the fast metals sig- 
nificantly contribute neither to the production of the X-ray lines 
from the ambient atoms nor to the bremsstrahlung continuum 
radiation (see Tatischeff et al. 1998 ). We further assume that the 
accelerated protons and a particles are in solar proportion, that 
is, C a /Cp = ane (see Eq. dA.7b ). 

In the calculations of the equilibrium spectra (N p and N a ), 
we neglect the nuclear destruction and catastrophic energy loss 
(e.g. interaction involving pion production) of the fast ions in 
the cloud. Indeed these processes are not important in compar- 
ison with the ionization losses below ~300 MeV nucleon" 1 ki- 
netic energy (see, e.g., Schlickeiser 2002) and most of the X-ray 
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Fig. C.l. Calculated differential equilibrium numbers of fast par- 
ticles (Nj\ solid lines) for the differential injection rate of pri- 
mary protons given by s-2 and E m \ n = 100 keV (dQ p /dt; 
dotted lines). Also shown are the differential production rates 
of secondary knock-on electrons (dQ e<s /dt; dashed lines), a): 
A = 10 21 cm" 2 ; b): A = 10 24 cm" 2 . The H density in the non- 
thermal X-ray production region is «h = 10 4 cm" 3 . The calcula- 
tions are normalized to a total power of 1 erg s" 1 injected by the 
primary LECR protons in this region. 



emission below 10 keV, which is the prime focus of the present 
work, is produced by ions in this low energy range (see Fig. [5j?). 



C. 1 . Secondary electron production 

We calculate the production of secondary electrons associated to 
the ionization of ambient H2 molecules and He atoms. The cor- 
responding differential ionization cross sections are obtained as 
in Tatischeff et al. ( fl998b from the work of Chu et al. ( fTMTT ). 
We neglect the production of secondary electrons and positrons 
that follows the production of charged pions in hadronic col- 
lisions. In fact, the corresponding electron and positron source 
functions can dominate the one of knock-on electrons only at 
energies > 10 MeV (Schlickeiser 2002), and these high-energy 
leptons are not important for the production of X-rays < 10 keV 
(seeFig.|5ji). 

Differential production rates of knock-on electrons are 
shown in Fig. IC.ll together with the corresponding equilibrium 
spectra of primary protons and secondary electrons. This figure 
illustrates the effects of changing the CR escape path length from 
A = 10 21 cm" 2 (panel a) to 10 24 cm -2 (panel b). In the first 
case, protons of energies up to 1 .4 MeV are stopped in the cloud, 
whereas in the second case the transition energy between proton 
stopping and escape is at 71 MeV. We see that above this transi- 
tion energy the equilibrium spectrum has a similar slope than the 
source spectrum, whereas at lower energies the equilibrium pro- 
ton distribution is harder due to the ionization losses. We can an- 
ticipate that the total X-ray production rate will be much higher 
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for the case A = 10 24 cm 2 , as a result of the higher proton 
number at equilibrium above a few MeV. 

C.2. X-ray continuum emission 

The X-ray continuum emission is due to inverse bremsstrahlung 
from the fast ions (the radiation of a single photon in the colli- 
sion of a high-speed ion with an electron effectively at rest) and 
classical bremsstrahlung from the secondary knock-on electrons. 
In the nonrelativistic domain, the bremsstrahlung produced by a 
proton of kinetic energy £ in a collision with a H atom at rest 
has the same cross section as that of an electron of kinetic energy 
(m e /m p )E in a collision with a stationary proton (m e and m p are 
the electron and proton masses, respectively). We calculate this 
cross section as in Sect. 3.2, but without taking the screening of 
the H nucleus by the bound electron into account. The cross sec- 
tion for interaction of a proton with a H atom is then multiplied 
by ( 1 +2aHe) to take the ambient He into account. For a particles, 
we replace the proton energy E by the energy per nucleon of the 
projectile and multiply the proton cross section by 4 to account 
for the nuclear charge dependence of the bremsstrahlung cross 
section. 

In the relativistic case, the cross section for proton inverse 
bremsstrahlung is different from the one for classical electron 
bremsstrahlung, owing to the appearance of angular and energy 
abberations in the transformation between the two rest frames 
of the interacting particles ( I laug 2003). We checked that these 
effects can be neglected in good approximation in the present 
work. 

In Fig. IC.2l we show two X-ray spectra corresponding to the 
particle equilibrium spectra presented in Fig. IC.ll We see that 
the continuum emission is dominated by inverse bremsstrahlung, 
which is a general rule independent of the model parameters (see 
Tatischeff et al. 119981 . We also see that, as expected, the X-ray 
production rate is much higher for A = 10 24 cm~ 2 than for A = 
10 21 cm" 2 , the difference being a factor of 22, 337 and 1054 at 
1,10, and 100 keV, respectively. 

C.3. X-ray line emission 

For producing X-ray lines from the ambient atoms, we take 
both the contribution from secondary electrons (see Eq. ( lB.3l l) 
and that from primary ions into account. The cross sections 
for K-shell ionization by proton and a-particle impacts are ex- 
tracted from the data library implemented by Pia et al. (2009) 
in the Geant4 toolkit for the simulation of particle induced X- 
ray emission (PIXE). We use the cross sections calculated in the 
ECPSSR theory with high-velocity corrections (Lapicki et al. 
2008). These cross sections are more accurate for mildly rela- 
tivistic projectiles than those previously employed by Tatischeff 
et al fllg98) . 

Proton and a-particle collisions with target atoms do not lead 
to significant line broadening effects caused by multiple simul- 
taneous ionizations. We thus adopt as before a width of 10 eV 
for all the lines (see Sect. 3.3). We note, however, that the X- 
ray lines produced by collisions of ions heavier than 4 He can be 
shifted by several tens of e V, significantly broadened and split up 
into several components (Garcia et al. 119731 1. For example, the 
Fe K« line produced by impacts of O ions of 1 .9 MeV nucleon 1 
is blueshifted by ~50 eV in comparison with the one produced 
by 5-MeV proton impacts, and has a FWHM of~100 eV (see 
Garcia et al. [19731 figure 3.55). 
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Fig. C.2. Calculated X-ray emission produced by LECR protons 
and or-particles interacting in a gas cloud of solar metallicity, 
for the differential injection rate of primary protons shown in 
Fig. 9. The contribution of accelerated a-particles is included as 
explained in the text, assuming in particular the solar abundance 
C a ICp = 0.0964. a): A = 10 21 crrT 2 ; b): A = 10 24 cm 4 . IB: in- 
verse bremsstrahlung; SEB: secondary electron bremsstrahlung. 
Photoelectric absorption is taken into account with a H column 
density of 10 22 cm -2 . 

The most intense line produced by LECR protons and a- 
particles is also the neutral Fe Ka line at 6.40 keV (Fig. IC.2) . 
This line has an EW of 2.31 and 0.80 keV in the spectra shown 
in Figs. IC.2b and b, respectively. The second strongest line in 
these spectra is the Si Ka line at 1 .74 keV; its EW is equal to 
309 eV in panel a and 152 eV in panel b. 
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